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Abstract 

This  research  develops  a  space-time  adaptive  processing  (STAP)  radar  model 
for  side-looking  (SL)  arrays  with  platform  maneuver  incorporation,  and  examines  ma¬ 
neuver  effects  on  Matched  Eilter  (ME)  performance  and  sample  support  homogeneity 
for  linear  and  planar  arrays  at  various  ranges.  A  primary  assumption  of  most  STAP 
models  is  the  absence  of  platform  maneuver  during  the  coherent  processing  interval 
(GPI).  This  research  abandons  the  platform  stationarity  assumption,  allowing  plat¬ 
form  pitch,  roll,  and  yaw  during  the  GPI.  Maneuver  rates  (set  to  800°/sec)  are  held 
constant  throughout  the  GPI  and  each  maneuver  type  is  examined  individually. 

All  maneuver  effects  have  a  negative  impact  on  both  ME  performance  and  sam¬ 
ple  support  homogeneity.  The  introduction  of  pitch  marginally  impacts  ME  Signal- 
to-Interface-plus-Noise-Ratio  (SINK)  due  to  the  lack  of  achieved  azimuth  resolution 
in  this  research,  but  decreases  available  homogenous  sample  support  24%  at  11  km 
and  20%  at  66  km.  Roll  impact  on  ME  performance  is  observed  only  in  the  planar 
array,  with  ME  SINR  dropping  4  dB  but  sample  support  homogeneity  is  unaffected. 
Yaw  effects  are  the  by  far  the  most  devastating,  dropping  output  SINR  15  dB  at  11 
km  and  66  km,  and  decreasing  available  homogeneous  sample  support  34%  at  11  km 
and  18%  at  66  km. 
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Space-Time  Adaptive  Processing 


For  Side-Looking  Arrays 
With  Platform  Maneuver 

I.  Introduction 

This  research’s  primary  objective  is  the  advancement  of  space-time  adaptive  pro¬ 
cessing  (STAP)  radar  data  models  for  the  side-looking  (SL)  array.  A  modeling 
process  often  requires  assumptions  that  limit  a  model’s  scope,  ideally  without  sacri¬ 
ficing  its  accuracy  or  integrity.  A  primary  assumption  of  most  STAP  models  is  the 
absence  of  platform  maneuver  during  the  coherent  processing  interval  (CPI).  Airborne 
radar  platforms  are  predominantly  modeled  as  stationary,  and  while  limited  published 
work  has  modeled  the  impact  of  that  assumption  on  the  forward-looking  (FL)  array, 
no  published  work  discovered  has  modeled  the  impact  on  the  SL  array.  This  research 
abandons  the  platform  stationarity  assumption,  allowing  platform  pitch,  roll,  and  yaw 
during  the  CPI.  The  primary  focus  of  this  research  is  the  maneuver  model  develop¬ 
ment,  based  on  the  popular  stationary  model  of  [3, 10].  Secondary  objectives  include 
characterizing  platform  maneuver  effects  on  the  clutter  environment,  examining  re¬ 
sulting  Matched  Filter  (MF)  output  Signal-to-Interference-plus-Noise- Ratio  (SINK), 
and  examining  maneuver  impact  on  data  homogeneity  and  sample  support. 

1 . 1  Motivation 

The  primary  motivation  for  this  research  is  contributing  an  additional  layer  of 
capability  and  robustness  to  the  popular  model  of  [3, 10].  Limited  published  work 
investigates  the  impact  of  platform  maneuver.  And  while  two  authors  have  examined 
the  FL  scenario,  the  clutter  models  used  were  their  own  and  not  widely  available  nor 
used.  This  research  offers  a  complete  framework  to  model  platform  maneuver  within 
the  context  of  the  popular  [3, 10]  model.  Future  researchers  may  use  this  framework 
to  examine  maneuver  effects  on  partially- adaptive  processing  techniques.  Others  may 
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use  it  to  explore  optimum  processing  techniques  less  susceptible  to  platform  maneuver. 
Others  may  convert  this  research  to  a  forward-looking  (FL)  array,  the  orientation  of 
most  fighter  aircraft  today.  This  thesis  lays  the  groundwork  for  future  investigations 
of  platform  maneuver.  Today’s  fiy-by-wire  aircraft  constantly  monitor  their  platform 
attitude,  offering  an  additional  potential  application  of  the  framework  this  research 
builds.  This  research  may  offer  potential  uses  in  the  space-based  radar  area  as  well, 
as  platform  maneuver  and  space-based  radar  are  inseparable. 

1 . 2  Organization 

The  stationary  model  of  [3, 10]  serves  as  the  foundation  for  this  research.  It  and 
STAP  fundamentals  are  discussed  in  Chapter  II.  The  maneuver  model  is  developed  in 
Chapter  III,  abandoning  the  stationary  platform  assumption  of  [3,10].  The  maneuver 
model  is  implemented  in  Chapter  IV,  and  maneuver  effects  on  the  clutter  environment, 
target,  and  radar  performance  are  examined.  Research  results,  improvements,  and 
future  relevant  research  incorporating  platform  maneuver  are  offered  in  Chapter  V. 

1.3  A  ssumptions 

Several  assumptions  within  the  stationary  model  of  [3,10]  remain  in  the  maneu¬ 
ver  model  developed.  The  waveform  is  considered  narrowband,  and  signals  are  co¬ 
phased  with  phase  shifts  and  not  time  delays.  The  waveform  is  assumed  Doppler  tol¬ 
erant,  making  the  time-frequency-autocorrelation-function  (TFACF)  unity.  Targets 
are  modeled  as  point  scatterers  with  constant  radar  cross  section  (RCS),  independent 
of  incidence  angle.  Maneuver  is  modeled  as  rotation  about  the  array  reference  ele¬ 
ment,  without  regard  to  array  location  on  the  platform.  Range  walk  effects  have  been 
previously  characterized  [5],  are  not  this  research’s  focus,  and  are  ignored.  Ground 
clutter  is  assumed  overwhelmingly  stronger  than  atmospheric  and  all  other  attenua¬ 
tion  effects,  and  is  the  only  factor  modeled.  Finally,  the  return  signal  is  assumed  a 
perfect  plane  wave,  with  no  phase  error  across  the  antenna  face. 
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II.  SL  Stationary  Model,  STAP  Fundamentals,  and 

Literature  Review 


Relevant  background  knowledge  consists  of  fonr  major  areas  discussed  in  this 
chapter:  an  airborne  radar  overview,  the  physical  airborne  radar  model  for 
a  stationary  platform,  space-time  adaptive  processing  (STAP)  concepts,  techniques, 
and  metrics,  and  previously  published  work  investigating  platform  maneuver  effects 
on  STAP  performance  for  the  forward-looking  (FL)  array. 

Concepts  and  details  of  the  popnlar  STAP  model,  developed  for  side-looking 
(SL)  linear  arrays  in  [10]  and  extended  to  planar  arrays  in  [3],  are  highlighted  in 
Section  2.2.  The  stationary  model  presented  serves  as  the  foundation  of  the  maneuver 
model  in  Chapter  III.  Major  components  of  the  model  include  the  geometry,  transmit 
signal,  receive  signal,  formatting,  antenna  pattern,  target  model,  noise  model,  and 
clutter  model.  A  major  assumption  of  this  popular  model,  abandoned  in  Chapter  III 
and  the  focus  of  this  thesis,  is  the  absence  of  platform  maneuver  during  the  coherent 
processing  interval  (CPI). 

STAP  concepts,  techniqnes,  and  metrics  are  highlighted  in  Section  2.3.  Included 
are  STAP  fundamentals,  fully  and  partially  adaptive  processing  techniques,  covariance 
estimation,  sample  support  homogeneity  issues,  and  performance  metrics. 

Finally,  previously  published  work  modeling  platform  maneuver  effects  on  FL 
STAP  are  discussed  in  Section  2.4.  Limited  published  work  for  the  SL  array  with 
platform  maneuver  is  available.  Two  significant  contribntors  for  the  FL  array  are 
Peter  G.  Richardson  [7,8]  and  George  M.  Herbert  [4]. 

Before  delving  into  the  stationary  model  development  and  STAP  concepts,  Sec¬ 
tion  2.1  offers  a  brief  overview  of  the  airborne  radar  problem. 

2.1  Airborne  Radar  Overview 

The  most  basic  and  fundamental  goal  of  a  radar  system  is  target  detection. 
Detection  requires  the  radar  to  distinguish  targets  from  the  interference  environment. 
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Pulsed  Doppler  radars  take  advantage  of  a  target’s  Doppler  frequency,  a  shift  due 
to  signal  time  compression  (or  dilation)  dependent  upon  the  target-radar  relative 
velocity  [6].  The  Doppler  frequency  fd  of  a  target  is  given  by  [9] 


Ao 


(2.1) 


where  Vr  is  the  relative  velocity  between  the  radar  and  target,  and  is  the  transmitted 
signal’s  wavelength.  In  the  ground-based  radar  scenario,  the  relative  velocity  between 
the  radar  and  ground  clutter  is  zero  since  neither  are  moving.  Therefore,  moving 
targets  that  produce  a  Doppler  shift  are  easily  separable  from  zero  Doppler  clutter. 
The  airborne  radar  does  not  possess  this  advantage.  Ground  clutter  returns  possess  a 
Doppler  frequency  in  the  airborne  radar  scenario  due  to  the  aircraft’s  velocity,  making 
the  task  of  separating  slow  moving  targets  from  clutter  more  difficult. 

A  basic  pulsed  Doppler  radar  system  transmits  a  series  of  pulses,  samples  the 
environment  at  specific  times  corresponding  to  range,  receives  the  return  pulses,  then 
adds  (or  equivalently,  integrates)  the  pulses.  The  method  of  integration  can  have  dra¬ 
matic  effects  on  detection  performance.  Non-adaptive  (NA)  techniques,  among  the 
earfiest  methods  of  processing  techniques,  simply  add  all  available  returns  without  any 
consideration  of  their  contents.  The  problem  arises  when  the  returns  contain  both 
the  desired  target  and  undesired  interference.  The  concept  of  adaptive-processing 
is  to  intelligently  weight  available  returns  to  minimize  interference  effects  and  max¬ 
imize  target  returns,  Signaf-to-Interference-pfus-Noise  Ratio  (SINK),  and  detection 
capabifity. 

In  order  to  intelligently  choose  these  adaptive  weights,  the  interference  envi¬ 
ronment  must  be  known  or  estimated.  The  airborne  radar  model  described  next  in 
Section  2.2  models  this  interference  environment  for  a  SL  antenna  array.  Section  2.3 
then  describes  fundamental  STAR  concepts,  techniques,  and  metrics  used  to  enhance 
target  detection. 
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Figure  2.1:  Antenna  array  geometry  for  stationary  model.  The  antenna  array  is 

oriented  in  the  x-z  plane.  Its  reference  element,  defined  as  the  npper-left  most  element, 
serves  as  the  origin  of  the  cartesian  and  radar  coordinate  systems.  The  elements  are 
uniformly  spaced  with  inter-element  spacings  of  dx  in  the  x  direction  and  dz  in  the 
^  direction.  The  N  azimnth  channels  are  numbered  0  to  —  1,  and  the  P  elevation 
channels  0  to  P  —  1. 

2.2  Airborne  Radar  Model  for  Stationary  Platform 

The  goal  of  the  radar  model  is  accurate  characterization  of  the  target  and  the 
interference  environment.  The  first  model  component  is  the  scenario  geometry.  Once 
the  geometry  is  defined,  the  transmit  and  receive  signals  are  described.  The  returns 
are  then  formatted  in  a  linear  algebra  framework.  Finally,  the  interference  (noise  and 
clntter)  environment  is  modeled.  The  following  sections  examine  each  of  these  pieces 
from  the  model  of  [3, 10]  in  detail. 

2.2.1  Geometry.  Consider  the  antenna  array  depicted  in  Figure  2.1.  The 
antenna  is  orientated  in  the  x-z  plane  with  N  horizontal  elements  and  P  elevation 
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elements.  The  horizontal  elements  are  numbered  from  zero  to  N  —  1  with  uniform 
inter-element  spacing  dx,  while  the  vertical  elements  are  numbered  zero  to  P  —  1  with 
uniform  inter-element  spacing  dz.  The  antenna  is  mounted  in  a  SL  configuration 
on  an  aircraft  above  the  ground,  as  depicted  in  Figure  2.2.  Boresight  is  defined  in 
the  y  direction.  The  aircraft  travels  at  speed  Va  in  the  x  direction  perpendicular 
to  the  y-z  plane.  A  target  is  located  at  an  elevation  angle  6  and  azimuth  angle  0. 
The  radar  elevation  angle  9  is  negative  when  measured  from  boresight  towards  the 
ground  and  the  radar  azimuth  angle  (p  is  positive  when  measured  from  boresight  to 
the  positive  x-axis  (toward  the  nose  of  the  plane).  The  upper  and  left  most  element 
is  the  first  element  to  receive  returns  from  a  target  at  positive  9  and  positive  0, 
and  is  therefore  designated  the  reference  element.  The  reference  element  defines  the 
origin  (0,0,0)  of  the  reference  cartesian  coordinate  system.  The  reference  element  and 
element  separations  dx  and  dz  define  individual  array  element  locations.  The  unit 
vector  from  the  reference  element  to  the  element  is  [3] 

dnp  =  —ndx'x.  —  pdzZ.  (2.2) 

The  vector  k  describes  the  unit  vector  to  a  target  at  (0,0)  and  a  transformation  from 
the  radar  coordinate  system  to  the  cartesian  coordinate  system.  Unit  vector  k  is  the 
projection  onto  the  cartesian  unit  vectors  x,  y,  and  z  [3], 

k  =  cos  9  sin  0x  -|-  cos  9  cos  0y  -|-  sin  0z.  (2.3) 

Equations  (2.2)  and  (2.3)  fully  describe  the  airborne  radar  SL  array  physical  picture. 
The  radar  transmit  and  receive  signals  are  next  characterized. 
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Reference  Element 


Figure  2.2:  Radar  platform  geometry.  The  array  is  mounted  on  an  aircraft  in  a  SL 
configuration  with  boresight  defined  in  the  y  direction.  The  aircraft  travels  at  speed 
Va  in  the  x  direction.  The  vector  k  is  the  unit  vector  from  the  reference  element  to 
the  target  at  (0,0). 

2.2.2  Transmit  Signal.  Construction  of  the  transmit  signal  model  begins 
with  a  single  pulse  (the  standard  envelope  function)  Up(t)  [3], 

{1  0  <  t  <  r 

0  otherwise , 

where  r  is  the  uncompressed  pulse  width.  Multiple  copies  of  Up{t)  shifted  by  integer 
multiples  of  the  pulse  repetition  interval  (PRI)  Tr  create  the  pulse  train  u{t), 


M-l 

u(t)  =  y  (t  —  mTr).  (2-4) 

m=0 
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The  summation  bounds  determine  the  CPI  length  of  M  pulses  or  equivalently  MT^ 
seconds.  Incorporation  of  a  signal  amplitude  Uj,  carrier  frequency  Wq,  and  random 
starting  phase  Lp  yields  the  complete  transmitted  pulsed  radar  signal  s{t), 

s{t)  =  (2.5) 

2.2.3  Received  Signal.  Upon  transmission  the  signal  propagates  through  free 
space  as  a  plane  wave,  reflects  off  a  target,  and  is  received  by  each  individual  array 
element  as  a  plane  wave.  Each  individual  array  element  contains  its  own  receiver 
chain.  The  received  signal  at  each  element  contains  three  major  differences  than 
the  original  transmit  signal:  a  Doppler  shift  fd  (assumed  equal  for  all  elements)  due 
to  relative  target  velocity  vp,  a  time  delay  due  to  the  round  trip  time  Tr  and 
differential  time  delay  and  a  received  amplitude  a^.  Incorporating  these  effects 
yields  the  received  signal  at  element  np  [3], 

Snp{t)  =  (2.6) 

The  receive  amplitude  includes  range  and  radar  cross  section  (RCS)  attenuation 
effects,  fully  discussed  in  [9].  The  time  delay  Tnp  is  the  sum  of  the  round  trip  time 
from  the  array  face  to  the  target,  Tr,  and  the  differential  time  from  the  reference 
element  to  each  individual  np*^  element, 

Tnp  =  Tr  +  T^p.  (2.7) 

A  narrowband  assumption  requires  a  signal  time  duration  greater  than  the  time  re¬ 
quired  to  propagate  across  the  face  of  the  array.  The  pulsewidth  must  therefore  meet 
the  requirement 

(2,8) 

c 
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where  \J {Nd^Y  +  {PdzY  represents  the  greatest  distance  across  the  array  face  and  c 
the  propagation  speed.  Under  the  narrowband  assumption  each  channel  receives  the 
return  signal  simultaneously  and  is  ignored  with  respect  to  time,  but  not  phase. 
Substituting  (2.7)  into  (2.6),  ignoring  the  differential  delay  in  the  pulse  train,  and 
rearranging  the  result  yields 


aj-u(t  Tfij 


gi27r(/o+/d)tg-i27r(/o+/d)Ti{g-j27r(/o+/dXpgjV 


(2.9) 


Complete  characterization  of  the  narrowband  receive  signal  therefore  requires  quan¬ 
tifying  the  receive  amplitude  and  the  signal  phase.  The  differential  delay  has 
units  of  time  and  represents  the  time  difference  of  signal  arrival  between  elements. 
The  corresponding  phase  difference  between  elements  must  be  found.  The  physical 
interpretation  of  is  the  element  to  element  physical  separation  in  the  k  direction 
divided  by  the  propagation  speed  c.  The  projection  of  the  unit  vector  to  the  target 
(k)  onto  the  physical  location  vector  of  the  elements  (d„p)  yields  that  element  to 
element  physical  separation,  resulting  in  the  expression 


np 


k(6>,  4>)  ■  dnp 
c 

—ndx  cos  0  sin  (/)  —  pdz  sin  9 
c 


(2.10) 

(2.11) 


Examining  the  third  exponential  term  in  Equation  2.9,  the  e  approximated 

as  e°  =  1,  since  0  for  any  expected  Doppler  values.  Multiplying  the  time 

difference  by  Ug  =  ‘27: fo  yields  the  element  to  element  phase  difference. 


/  o  £  ^  +  pdz  sin  9 

^o^rip  2::  j  Q 


(2.12) 
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A  spatial  frequency  {)  defines  the  return  signal  element  to  element  phase  behavior, 


k(6*,  0)  •  d 


np 


Xo 

dx  cos  0  sin  0  ^  dz  sin  6 


\n 


\n 


The  x-axis  and  z-axis  spatial  frequencies  are  given  by 

dx  cos  9  sin  (p 


and 

dz  sin  9 

^ - , 

respectively.  Equation  (2.12)  is  then  rewritten  as 

=  -27T{ni9x  +  pidz). 

Incorporating  Equation  (2.17)  into  Equation  (2.6)  results  in 

Snpit)  =  aru{t  -  TR)e^'2^(A+A)ie-i2^(/o+/dTRgi/ci<pgi2^M.+ptfdeiv^ 


(2.13) 

(2.14) 


(2.16) 


(2.16) 


(2.17) 


(2.18) 


Since  /<,,  fd,  and  Tr  are  not  functions  of  time  and  are  constants  for  a  given  Snp(t),  the 
^-j2TT{fo+fd)TR  is  rolled  into  the  random  phase  ip.  The  complete  receive  signal 

therefore  becomes  [3] 


Snp{t)  =  aru{t  -  (2.19) 


Characterization  of  the  signal  at  each  receiver  channel  is  now  complete.  Each  receiver 
channel  performs  a  matched  filtering  operation  and  analog  to  digital  conversion  after 
stripping  the  carrier  frequency,  as  depicted  in  Eigure  2.3.  The  resulting  signal  is  [3] 


X 


mnp 


(2.20) 
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s{t) 


Figure  2.3:  Element  channel  block  diagram.  The  received  signal  s(t)  passes  through 
a  low  noise  ampliher,  is  stripped  of  its  carrier  frequency,  match  hltered,  then  converted 
to  the  digital  signal  x^np- 


where  at  is  a  complex  amplitude  accounting  for  amplitude  and  phase  (f>,  m  is  the 
pulse  within  the  CPI  corresponding  to  the  single  time  instance  (range  cell)  the  signal 
was  sampled,  and  ui  is  the  normalized  doppler  frequency  given  by 


(2.21) 


Complex  amplitude  at  accounts  for  all  range  effects  of  the  radar  range  equation  [9]. 
The  matched  filtering  step  glossed  over  in  Equation  (2.20)  assumed  a  Doppler  tolerant 
waveform.  The  step  is  valid  as  long  as  [5] 


Vr 


(2.22) 


where  /„  is  the  transmit  frequency  and  B  is  the  receiver  bandwidth.  An  unrealistic 
relative  velocity  on  the  order  of  90  x  10^  meters  per  second  would  be  required  to 
invalidate  the  step  for  the  parameters  investigated  in  this  research. 


2.2.4  Formatting.  For  a  given  range  cell,  the  radar  has  returns  from  M 
pulses,  N  azimuth  channels,  and  P  elevation  channels  for  a  total  of  MNP  returns. 
These  MNP  returns  are  arranged  in  a  compact  manner  suitable  for  linear  algebra  op¬ 
erations  to  integrate  the  returns.  Arrangement  is  accomplished  in  [3]  through  spatial 
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(azimuth  and  elevation)  and  temporal  steering  vectors.  The  azimuth  steering  vector 
describes  the  element  to  element  phase  progression  due  to  the  physical  separation 
among  N  azimuth  channels.  The  elevation  steering  vector  describes  the  element  to 
element  phase  progression  due  to  the  physical  separation  among  P  elevation  chan¬ 
nels.  And  the  temporal  steering  vector  describes  the  pulse  to  pulse  phase  progression 
due  to  the  time  separation  of  arrival  of  the  M  pulses.  The  three  steering  vectors  are 
arranged  to  represent  the  element  to  element,  pulse  to  pulse  phase  progression  and 
amplitudes  of  the  MNP  returns.  This  MNP  x  1  vector  is  the  steering  vector  Vj  to 
location  (d,0). 

The  spatial  steering  vectors,  a(dx)  and  e(dz),  describe  the  element  to  element 
phase  progression  of  the  return.  With  the  upper-left  most  element  labeled  the  0 
reference  element  and  the  last  of  N  azimuth  channels  labeled  —  1,  the  element  to 
element  phase  progression  in  azimuth  becomes, 

a(d^)  =  _  ^2.23) 

With  the  upper-left  most  element  labeled  the  0  reference  element  and  the  last  of  P 
elevation  channels  labeled  P—1,  the  element  to  element  phase  progression  in  elevation 
becomes, 

e(d^)  =  (2.24) 

The  temporal  steering  vector,  b((h),  describes  the  pulse  to  pulse  phase  progression  of 
a  return  at  a  single  normalized  doppler  frequency,  u.  The  pulse  to  pulse  phase  delay 
due  to  Doppler  is  a  linear  function  of  a;  if  a)  is  constant, 

b(cD)  =  ^2.25) 
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The  spatial  and  temporal  steering  vectors  are  combined  into  the  single  steering  vector 
Vt  using  the  Kronecker  product  0  [3], 

vt(d,  0)  =  6(^2)  0  b(a))  0  a(d^).  (2.26) 

The  steering  vector  vj  captures  the  element  to  element,  pulse  to  pulse,  phase  progres¬ 
sion  of  a  signal  return  from  a  single  0,  single  (f>,  and  single  a).  Incorporating  the  signal 
amplitude  at  yields  the  element  to  element,  pulse  to  pulse,  return  signal  (magnitude 
and  phase)  called  the  space-time  snapshot  and  given  by. 


Xt  =  (2.27) 

Actual  airborne  radar  returns  may  contain  multiple  targets,  clutter,  noise,  jam¬ 
ming,  and  other  undesirable  effects.  The  complete  space-time  snapshot  therefore 
contains  each  of  these  factors.  However,  the  only  effects  considered  in  this  thesis 
are  target,  noise,  and  clutter.  The  complete  snapshot  therefore  contains  the  target 
snapshot  Xt^  noise  snapshot  and  the  clutter  snapshot  Xc-  The  target  present 
scenario  is  called  the  Hi  hypothesis, 

Xm  =  Xt  +  Xn  +  Xc-  (2.28) 

The  target  absent  scenario  is  called  the  Hq  hypothesis  and  contains  only  undesirable 
effects, 

Xho  =  Xn  +  Xc-  (2-29) 

The  signal  model  and  formatting  are  now  in  place.  The  overall  goal  remains 
characterization  of  the  target  and  the  interference  (thermal  noise  and  clutter  for  this 
research)  environment.  The  next  step  is  using  the  formatting  framework  developed 
to  accurately  characterize  target,  noise,  and  clutter  space-time  snapshots.  These 
snapshots  then  create  the  models  themselves.  The  antenna  pattern  is  an  important 
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component  of  the  target  and  clutter  models  as  it  directly  relates  to  the  return  ampli¬ 
tudes.  The  antenna  pattern  is  defined  next,  followed  by  the  target,  noise,  and  clutter 
model. 


2.2.5  Antenna  Pattern.  The  magnitude  of  target  and  clutter  returns  depend 
on  the  power  with  which  they  are  illuminated.  That  illumination  depends  on  the 
array’s  antenna  pattern.  Antenna  pattern  development  requires  an  element  pattern 
and  an  array  factor.  Multiplication  of  the  element  pattern  and  array  factor  yield  the 
complete  antenna  pattern. 

The  array  elements  are  assumed  identical  in  power  and  voltage  patterns.  A 
cosine  pattern  in  both  azimuth  and  elevation  with  a  backlobe  attenuation  factor  of  bg 
fully  defines  the  element  voltage  pattern  [3], 


cos  0  cos  0  —90°  <  6,(j)  <  90° 

beCosdcostp  90°  >  0,0  >  270° 


f  (»,<!>) 


Backlobe  attenuation  represents  a  screen  behind  the  array  [9].  The  array  factor  for 
an  array  steered  to  zero  azimuth  and  zero  elevation  is  [1], 


P-l  JV-l 


(2.30) 


p=0  n=0 


Manipulation  into  an  easily  implemented  form  is  accomplished  in  [3]  using  geometric 
series  properties, 


1T(0,0) 
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Undefined  points  in  the  form  given  by  Equation  (2.31)  are  determined  in  [3]  using 
L’Hopital’s  rule, 


W{9,Q)  = 


sin(z^Psind) 
sin(5;Jsind)  ’ 

sin  A/"  sin  0) 
gin  (2^  sin  0)  ’ 


(2.32) 

(2.33) 


and 


tU(0,0)  =  NP. 


(2.34) 


The  product  of  the  element  pattern  and  array  factor  yields  the  spatial  array  pattern 
in  volts, 

P(0,0)  =  lU(0,0)m0),  (2.35) 

or  can  be  expressed  in  watts, 


G'(0,0)  =  \w{eA)f{0A)? 

=  |lU(d,(/))|2p(0,0). 


(2.36) 

(2.37) 


Equation  (2.37)  describes  the  radiated  power  in  any  {9^<p)  direction,  represent¬ 
ing  the  final  component  needed  to  accurately  capture  a  target  or  clutter  space-time 
snapshot.  The  target,  noise,  and  clutter  models  are  now  ready  to  be  defined. 

2.2.6  Target  Model.  The  MNP  x  1  target  snapshot  describes  the  target 
return  amplitude,  aj,  and  pulse  to  pulse,  element  to  element  phase  progression  for  a 
target  at  a  single  d,  single  0,  and  single  7), 

Xt  =  (2.38) 

=  ate{6z)  ®  b(a;)  0  a('da,).  (2.39) 

Equations  (2.23)  and  (2.24)  define  the  spatial  steering  vectors  ail'd z)  and  e(d^),  while 
Equations  (2.15)  and  (2.16)  define  the  spatial  frequencies  and Equation  (2.25) 
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defines  the  temporal  steering  vector  b(a;).  The  normalized  Doppler  frequency  uj  is 
defined  [9], 


(2.40) 


where  is  the  relative  velocity  between  the  aircraft  and  target,  A  is  the  radar  transmit 
wavelength,  and  fj.  is  the  pulse  repetition  frequency  (PRF).  The  relative  velocity 
between  two  objects  is  the  time  rate  of  change  of  the  distance  between  them, 


Vy  = 


dR 

dt 


(2.41) 


This  time  rate  of  change  is  the  projection  of  the  aircraft  velocity  vector  onto  k 
plus  the  projection  of  the  target  velocity  vector  Vj  onto  k, 


Vr 


k  •  Va  +  k  ■  Vt 


(2.42) 


In  the  stationary  model  of  [3, 10],  aircraft  velocity  is  restricted  to  the  x  direction  and 
thus  given  by 

Va  =  Va±.  (2.43) 

Target  velocity  may  include  x,  y,  and  z  components,  denoted  uf,  and  respec¬ 
tively, 

vt  =  <x  +  ufy  +  v^z.  (2.44) 

The  normalized  Doppler  frequency  for  a  target  at  (0,0)  is  therefore 

{va  +  v^)  cos  0  sin  0  -|-  uf  cos  0  cos  0  +  Uj  sin  0 

- AA - 

Any  change  in  the  relative  velocity  requires  a  recalculation  of  the  normalized  Doppler 
frequency.  In  the  stationary  model  of  [3,10],  the  relative  velocity  and  subsequently 
Doppler  frequency  of  the  target  remains  constant  throughout  the  CPI.  Complete  char¬ 
acterization  of  Vt  is  therefore  not  required.  In  the  maneuver  model  of  Chapter  III, 


2.45 
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relative  velocity  and  subsequently  Doppler  frequency  are  not  constant  and  character¬ 
ization  of  Vj  is  required. 

The  remaining  component  of  the  target  space-time  snapshot  is  the  target  am¬ 
plitude,  at-  A  target  to  noise  ratio  is  dehned  [3], 

,  P,G(0,4>)g(0,4>)>fa,  . 

where  Pt  is  the  transmit  power,  G  is  the  array  gain  on  transmit,  g  is  the  element  gain 
upon  receive,  at  is  the  target  RCS,  R  is  the  range  to  target,  and  NoBLg  accounts  for 
system  noise  and  losses.  The  target  RCS  is  assumed  independent  of  aspect  angle  and 
constant  throughout  the  CPI.  The  complex  amplitude  at  obeys 

£{atat}  =  a'^^t:  (2.47) 

where  a^  is  the  receiver  thermal  noise  power  and  £{  }  is  the  expected  value  oper¬ 

ator.  The  thermal  noise  power  is  defined  by  the  receiver  operating  temperature  and 
bandwidth, 

CT^  =  kToB  =  NoB,  (2.48) 

where  k  is  Boltzmann’s  constant.  In  the  implementation  of  the  stationary  model 
[3, 10],  ^t  is  set  equal  to  one.  In  Chapter  Ill’s  maneuver  model,  ^t  is  recalculated  for 
every  change  in  platform  orientation. 

The  target  space-time  snapshot  is  now  described.  Equation  (2.46)  and  (2.47) 
produce  the  target  return  amplitude.  The  target  return  phase  progression  is  described 
by  Equations  (2.23),  (2.24)  and  (2.25).  The  remaining  components  of  the  stationary 
model  are  the  noise  and  clutter  characterizations. 

2.2.7  Noise  Model.  Only  internally  generated  receiver  noise  is  considered 
and  each  element  has  its  own  receiver.  The  noise  is  therefore  temporally  and  spatially 
white  and  uncorrelated.  The  noise  correlation  matrix  is  found  with  the  expected  value 
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of  the  outer  product  of  the  noise  snapshot  and  its  Hermitian, 


R..  =  (2.49) 

The  subscript  n  throughout  Equation  (2.49)  represents  “noise”,  not  to  be  confused 
with  array  elements.  The  lack  of  correlation  temporally  and  spatially  reduces  the  non¬ 
diagonal  values  of  S{XnXn}  zero.  The  diagonal  values  take  the  form  ^{XmnpXmnp}; 
where  the  subscript  n  here  represents  the  element  number.  The  notation  is  somewhat 
ambiguous,  but  follows  the  convention  of  [3, 10].  The  expected  value  of  the  thermal 
noise  voltage  return  times  its  conjugate  is  simply  the  thermal  noise  power,  cr^.  The 
resulting  noise  correlation  matrix  is  the  thermal  noise  power  times  a  MNP  x  MNP 
identity  matrix  [3], 

R-n  =  (2.50) 

The  presence  of  thermal  noise  is  vital  to  STAP  applications.  Thermal  noise  ensures 

the  interference  covariance  matrix  is  invertible.  The  final  piece  of  that  interference 
environment  for  this  research  is  the  clutter  model. 

2.2.8  Clutter  Model.  A  MNP  x  1  clutter  snapshot  describes  the  target 
return  amplitude,  ckc,  and  pulse  to  pulse,  element  to  element  phase  progression  for  a 
clutter  patch  at  a  single  0,  single  0,  and  single  (h, 

Ac  =  «cVc  (2.51) 

=  a;ce(d2)  ®  b((h)  ®  a(da;).  (2.52) 

Equations  (2.23)  and  (2.24)  define  the  spatial  steering  vectors  and  e{'dz),  while 
Equations  (2.15)  and  (2.16)  define  the  spatial  frequencies  dx  and  1?^.  Equation  (2.25) 
defines  the  temporal  steering  vector  b((x;).  The  normalized  Doppler  frequency  lD  is  de¬ 
fined  in  Equation  (2.21.  The  clutter  patch  is  stationary,  meaning  the  relative  velocity 
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depends  solely  on  the  aircraft  velocity, 


Vr  =  k  •  Va 

=  Va  COS  6  sin  (p 


(2.53) 

(2.54) 


The  normalized  Doppler  frequency  of  a  clutter  patch  at  {O,^)  is  therefore 


uj  —  2 


Va  COS  9  sin  p 
A/r 


(2.55) 


Next,  the  clutter  amplitude  «c  must  be  characterized.  The  elevation  angle,  grazing 
angle,  and  physical  attributes  of  the  clutter  patch  must  first  be  described.  Referring 
to  Figure  2.4,  the  spherical  Earth  model  uses  a  4/3  effective  radius  [3, 10], 


Qjq 


(2.56) 


where  Ve  represents  the  Earth’s  true  radius.  The  effective  radius  model  assumes  a 
uniform  vertical  gradient  of  atmospheric  refractive  index,  accounting  for  the  atmo¬ 
spheric  refraction  of  the  propagated  waveform.  The  elevation  angle  represents  the 
angle  between  the  line  to  a  clutter  patch  at  range  Rc  and  boresight,  and  is  found 
using  the  law  of  cosines  [3] , 


9^-  sin-^ 


/ +  ha{ha  2ae)A 
\  2Rc{ha  +  Cle)  ) 


(2.57) 


The  grazing  angle  represents  the  angle  between  the  line  tangential  to  the  surface  at 
Rc  and  the  line  from  the  reference  element  to  the  patch  at  Rc  [3], 


9ai9a  T  2(2g) 
2RcCic 


(2.58) 


The  horizon  range  is  defined  at  the  range  at  which  the  grazing  angle  is  zero. 


Rh  5/^1  T 


(2.59) 
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Figure  2.4:  Clutter  model  Earth  geometry.  The  Earth  is  modeled  with  a  4/3 

effective  radius  The  aircraft  is  located  at  height  ha-  The  elevation  angle 

0  represents  the  angle  from  boresight  to  patches  at  range  Re-  The  grazing  angle  xp 
represents  the  angle  between  the  line  tangential  to  the  surface  at  Rc  and  the  line  from 
the  reference  element  to  the  patch  at  Rc- 

The  ground  plane,  located  in  the  x-y  plane  at  height  —ha  and  depicted  in  Eigure  2.5, 
is  divided  into  i  x  k  individual  clutter  patches.  The  dimensions  of  these  patches  are 
determined  by  the  radar  operating  parameters.  The  number  of  clutter  patches,  Nc, 
determines  the  clutter  patch  angular  extent,  A0, 

A(f)  =  — .  (2.60) 

The  radar  range  resolution,  ^  for  the  uncompressed  pulse,  determines  the  clutter 
patch  range  extent,  AR, 

(2.61) 
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Figure  2.5:  Clutter  ground  plane.  The  ground  plane  is  divided  into  ixk  individual 
clutter  patches  of  angular  extent  A0  =  2tt/Nc  and  range  extent  AR  =  crj^.  The 
radar’s  unambiguous  range  Ru  and  horizon  range  Rh  determine  the  range  ambiguities 
Nr-  The  range  cell  under  test  at  Rc  is  the  unambiguous  ring  and  therefore  labeled 
z  =  0.  Ambiguous  rings  are  located  at  integer  multiples  of  the  unambiguous  ring.  In 
the  condensed  illustration  above,  the  unambiguous  ring  at  Rc  is  the  second  range  ring. 
The  first  ambiguous  ring,  labeled  z  =  1,  is  therefore  the  fourth  range  ring.  The  next 
ambiguous  ring  (not  shown),  is  the  sixth  range  ring.  Ambiguous  rings  are  counted 
until  the  horizon  range  Rh  is  reached. 

The  number  of  range  ambiguities,  are  given  by, 

RcNr  >  Rh 
RcNr  <  Rh  • 
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With  the  angles  and  physical  attributes  defined,  the  clutter  patch  amplitude  is 
quantifiable.  The  Clutter-to-Noise- Ratio  (CNR),  is  defined  [3], 

_  P,G{0,  me.  <!>)>? 

(4,7r)>NBLA(0)i  '  ^  ’ 

where  Pt,  G,  g,  A,  NBLg,  and  Rc  were  previously  defined.  Unlike  the  target  model, 
the  clutter  patch  RCS,  (Tc,  varies  with  respect  to  the  angle  of  incident  illumination. 
The  clutter  patch  RCS  is  a  function  of  the  clutter  patch  area  ARA(l)R,  reflectivity 
j  sin 'ip  [9],  range  Rc,  and  grazing  angle  'ip, 

(Tc  =  7  sin  'ipARAipR  sec  'ip.  (2.63) 

The  patch  reflectivity  is  modeled  as  constant  gamma  (7)  [9].  The  complex  amplitude 
Oc  of  a  smgle  clutter  patch  obeys  [3] 

£{acal}  =  (T^^c.  (2.64) 

But  the  antenna  array  radiates  in  all  directions,  illuminating  every  ixk  clutter  patch. 
Therefore,  the  complete  clutter  space-time  snapshot  is  the  sum  of  i  x  k  clutter  patch 
snapshots  [3], 

Nr-l  Nc-1 

z{i)^  ®  b(^i/c)  ®  ^(^x(i,k))  ■  (2.65) 

i=0  k=0 

And  the  complex  amplitude  aik  obeys  [3] 

£{aika*i]  =  a^CikSi-j5k-i,  (2.66) 
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meaning  the  returns  from  different  clutter  patches  are  uncorrelated.  The  final  piece 
of  the  clutter  model  is  the  clutter  covariance  matrix  from  [3] , 

R'c  =  ^{XcXc) 

Nr-1  Nc-1 

=  cr^  O  b((h)b((h)^  O  a(d^)a(d^)^.  (2.67) 

i=0  k=0 

Environment  characterization  is  now  complete.  The  target,  clutter,  and  noise 
models  are  described  and  the  stationary  SL  radar  is  able  to  be  simulated.  The  frame¬ 
work  and  concepts  of  this  stationary  model  serve  as  the  foundation  of  the  maneuver 
model  in  Chapter  III.  There  are  two  final  areas  to  visit  before  maneuver  model  de¬ 
velopment  begins.  The  first  covers  an  overview  of  STAP  fundamentals,  followed  by  a 
succinct  review  of  previously  published  work  modeling  platform  maneuver  effects  on 
STAP. 


2.3  Fundamental  STAP  Concepts  and  Metrics 

This  section  discusses  four  fundamental  areas  critical  to  understanding  the  ob¬ 
jectives  of  this  research.  The  first  is  general  STAP  concepts,  discussed  in  Section 
2.3.1.  The  second  is  sample  support  homogeneity  discussed  in  Section  2.3.2.  The 
third  is  clutter  spectrum  estimation  techniques  discussed  in  Section  2.3.3.  And  the 
final  area  is  STAP  metrics  discussed  in  Section  2.3.4. 

2.3.1  General  STAP  Coneepts  .  The  motivation  for  STAP,  and  radar  de¬ 
tection  in  general,  is  maximizing  target  detection  and  minimizing  interference  effects. 
The  general  definition  for  a  radar  filtered  output  is 

y  =  w^x,  (2.68) 

where  %  is  all  available  returns,  and  w  is  the  weight  vector  describing  how  those 
returns  are  combined.  The  stationary  model  of  [3,10]  reviewed  in  Section  2.2  provides 
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a  method  to  model  those  radar  returns  x-  The  scalar  output  y  is  then  threshold 
detected  to  determine  if  a  target  is  present.  The  choice  of  weight  vectors  w  has 
a  signihcant  impact  on  detection  performance.  A  poor  selection  of  weights  fails  to 
mitigate  the  interference  environment  and  losses  in  SINK  and  detection  are  incurred. 
Non- adaptive  processing  refers  to  the  weight  vector  being  the  steering  vector  v  to  the 
radar  look  location  (6*,0),  given  by  Equation  (2.26),  i.e. 

l/na  =  (2.69) 

STAP  seeks  to  optimize  SINK  and  target  detection  by  intelligently  selecting  the  weight 
vector  w  based  on  the  interference  statistics  of  the  environment.  The  model  of  [3, 10] 
characterizes  the  interference  by, 


R,  —  R-n  T  Rc)  (2.70) 

where  R„  and  Rc  are  defined  by  Equations  (2.49)  and  (2.67)  respectively.  The  opti¬ 
mum  weight  vector  is  known  as  the  Matched  Eilter  (ME)  and  given  by, 

Wmf  =  R“^V,  (2.71) 

where  R  is  the  known  interference  covariance  and  v  is  the  steering  vector  to  the 
target  location.  Unfortunately,  the  interference  environment  is  never  known  and 
must  be  estimated.  The  interference  estimate  is  accomplished  using  snapshots  at 
range  cells  (referred  to  as  “sample  support  vectors”)  surrounding  the  Range-cell- 
Under-Test  (RUT), 

(2.72) 

i=l 

2.3.2  Sample  Support  Homogeneity  .  The  accuracy  of  the  interference  esti¬ 
mate  depends  on  the  sample  support  vectors  used.  An  underlying  premise  of  inter- 
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ference  estimation  is  independent-identically-distributed  (i.i.d.)  sample  support  data. 
I.i.d.  sample  support  means  all  support  vectors  have  the  same  probability  distribu¬ 
tion  and  are  mutually  independent.  Non-homogeneities  within  the  sample  snpport 
vectors  cause  a  poor  interference  estimate.  A  Non-Homogeneity  Detector  (NHD)  is  a 
method  to  select  the  most  homogenous  sample  support  vectors,  thereby  minimizing 
degradation  cansed  by  non-homogeneons  support  data.  For  this  research,  NHD  will 
be  used  to  evaluate  platform  maneuver  induced  heterogeneities  in  the  support  data. 
The  generalized  inner-product. 


GIPk  —  X^H-RUTXk-  (2.73) 

provides  a  test  for  measuring  the  homogeneity  of  sample  support  data  [11].  Using 
the  known  covariance  of  the  RUT,  the  GIF  compares  the  magnitude  and  phase  of 
snapshots  at  snrrounding  range  cells  (designated  k).  The  larger  the  GIF  variance  and 
range  across  a  set  of  sample  support  vectors,  the  greater  the  non- homogeneity. 

One  of  the  reasons  NHD  is  important  is  the  unavailability  of  sufficient  sample 
snpport  for  interference  estimation.  Real  world  radars  have  MNP  degrees  of  freedom 
(DOF)  on  the  order  of  several  thousand.  Reed’s  rule  [2]  states  the  number  of  sam¬ 
ple  support  range  cells  required  to  achieve  SINR  within  3  dB  of  MF  performance  is 
two  times  the  DOF.  Rarely  are  there  anywhere  near  2MNP  range  cells  available  for 
fully  adaptive  processing.  The  lack  of  snfficient  sample  snpport,  as  well  as  prohibitive 
processing  reqnirements  for  fully  adaptive  processing,  necessitates  partially  adaptive 
processing  techniques  like  Factored-Time  Space  (FTS)  and  Joint-Domain-Localized 
(JDL).  While  these  techniques  are  not  examined  in  this  research,  the  manenver  in- 
dnced  heterogeneities  modeled  are  important  to  these  partially  adaptive  processes. 
FTS  requires  a  minimum  of  2N  i.i.d.  sample  support  vectors.  JDL  requires  2rjarjbrje 
i.i.d.  sample  support  vectors,  where  rja,  rib,  and  rje  are  the  DOF  used.  Therefore, 
maneuver  effects  on  data  homogeneity  will  influence  iid  sample  snpport  availability 
and  selection. 
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The  development  of  the  maneuver  model  is  the  primary  research  focus,  fol¬ 
lowed  by  examination  of  maneuver  effects  on  MF  performance  and  data  homogeneity. 
While  the  GIF  provides  an  evaluation  of  sample  support  homogeneity,  clutter  spec¬ 
tral  estimation  allows  insight  into  the  interference  environment  and  potential  impact 
on  detection  performance.  Two  spectral  estimation  techniques  are  employed  in  this 
effort. 


2.3.3  Spectral  Estimation  .  The  hrst  spectral  technique  used  is  the  low- 
resolution  Signal  Match  (SM).  The  SM  technique  steers  through  the  interference  data 
with  the  steering  vector  v,  producing  a  maximum  when  the  steering  vector  and  inter¬ 
ference  data  match.  The  technique  produces  spurious  sidelobes  due  to  the  3-D  discrete 
Fourier  transform  equivalence  of  the  MNP  x  1  steering  vector  used  in  the  [3,10]  model 
and  this  research.  These  sidelobes  are  not  a  part  of  the  clutter  model,  but  are  what 
the  radar  “sees”  when  processing  interference  data  with  the  the  steering  vector  and 
are  thus  well  suited  for  detection  purposes.  The  SM  spectral  estimate  is  given  by  [5] 


v^Rv 

■ 


(2.74) 


The  second  technique  is  the  high- resolution  Minimum- Variance  Estimate  (MVE).  The 
MVE  technique,  given  by  [5] 


Emve  =  (V^R  V)  \  (2.75) 

and  seeks  minimums  and  provides  the  location  of  interference  in  the  environment. 
MVE  amplitude  information  is  useless  as  it  only  provides  interference  source  location. 
The  MVE  leads  to  realistic  spectral  estimates  if  the  signal  is  composed  of  single 
spectral  lines,  but  unrealistic  estimates  if  it  is  a  continuous  spectrum  [5] . 

2.3.4  Metrics  .  Output  SINK  is  one  of  the  most  common  metrics  for  radar 
performance.  As  its  name  suggests,  it  is  a  ratio  of  the  desired  target  power  over  the 
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undesired  interference  power.  Output  SINK  is  given  by 


SINR  = 


I  p 

w^^Rw  ’ 


(2.76) 


where  is  the  steering  vector  to  the  target,  is  the  Signal-to-Noise- Ratio,  and  w  is 
the  weight  vector  applied.  Another  useful  metric  is  SINR  Loss, 


SINRi  = 


SINR 

itMNP' 


(2.77) 


SINR  Loss  references  the  output  SINR  to  the  thermal  noise  only  situation.  SINR  Loss 
is  a  convenient  tool  to  compare  results  from  scenarios  with  different  MNP  values  since 
it  is  normalized  by  MNP.  The  final  metric  used  is  the  Improvement  Factor.  The 
Improvement  factor  is  the  ratio  of  the  output  SINR  to  input  SINR.  The  larger  the 
improvement  factor,  the  better  the  radar  is  adapting  to  the  interference  environment. 
The  improvement  factor  is  given  by. 


SINR^t 
1  +  ^ 


where  is  the  Clutter-to-Noise-Ratio. 

These  metrics  are  used  in  Chapter  IV  to  analyze  the  effects  of  platform  maneu¬ 
ver. 


2.4  Literature  Review 

Peter  G.  Richardson  first  explored  platform  maneuver  effects  on  STAP  in  1997  [7]. 
Richardson  focused  on  roll  effects  for  the  forward-looking  (FL)  circular  planar  array 
and  their  implications  to  airborne  intercept  applications.  The  aircraft  velocity  vector 
was  fixed  in  [7] ,  while  the  maneuver  model  developed  in  Chapter  III  has  the  ability  to 
model  a  fixed  or  moving  aircraft  velocity  vector.  Richardson  observed  an  11  dB  drop 
in  SINR  when  roll  was  introduced,  but  predicted  motion  compensated  steering  would 
reduce  the  degradation.  Motion  compensated  steering  is  not  investigated  in  this  re- 
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search,  since  the  main  objectives  are  model  development  and  homogeneous  sample 
support  investigation.  However,  the  maneuver  model  developed  has  the  ability  to 
fully  simulate  motion  compensated  steering. 

In  1998  and  1999,  Richardson  extended  his  original  work  for  the  FL  circular 
planar  array  and  determined  roll  rates  as  small  as  0.005  degrees  per  PRI  significantly 
impacted  slow  moving  target  detection  [8] .  These  results  were  for  slow  moving  ground 
targets  with  the  antenna  steered  to  the  ground  location.  He  also  showed  allowing 
aircraft  velocity  vector  movement  had  minimal  impact  on  the  results.  This  minimum 
impact  was  also  observed  in  the  Chapter  HI  maneuver  model  developed.  A  significant 
portion  of  Richardson’s  work  in  [7,8]  was  roll  effects  on  jammer  rejection  performance 
in  the  FL  scenario. 

Minimal  research  was  published  after  Richardson’s  until  George  M.  Herbert 
revisited  the  topic  in  2005  [4].  Herbert  investigated  roll,  pitch,  and  yaw  effects  for 
the  FL  planar  array.  Using  similar  simulation  parameters  to  Richardson  in  [7]  and 
[8],  Herbert  concluded  roll  effects  were  not  nearly  as  significant  as  [7,8]  suggested. 
Herbert  attributed  the  pessimistic  results  of  [7, 8]  to  discontinuities  in  the  antenna 
pattern  and  a  simplistic  clutter  model.  Additionally,  Herbert  investigated  motion 
compensated  steering  and  determined  it  indeed  compensated  the  negative  effects  of 
platform  maneuver  as  Richardson  suggested. 

The  research  to  follow  builds  a  framework  to  model  pitch,  roll,  and  yaw  platform 
maneuver  for  the  SL  planar  array.  Several  deviations  from  the  previously  published 
work  of  Richardson  and  Herbert  exist.  The  most  obvious  is  the  SL  array  is  investigated 
instead  of  the  FL.  Next,  the  STAP  model  and  framework  is  based  on  that  of  [10] 
and  [3].  Third,  aircraft  velocity  vector  movement  is  allowed.  Motion  compensated 
steering  will  not  be  explored  and  the  antenna  mainbeam  is  steered  to  boresight  in  this 
research.  And  finally,  both  linear  and  planar  arrays  are  modeled  instead  of  circular 
or  linear  arrays.  Rotation  rates,  total  degree  of  rotations,  and  CPI  lengths  are  similar 
to  [4,7,8]. 
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The  overall  goal  of  both  the  stationary  model  of  [3, 10]  and  the  maneuver  model 
developed  in  this  research  are  accurate  characterization  of  the  interference  environ¬ 
ment.  The  approaches  taken  are  similar,  but  the  paths  are  slightly  different.  Chapter 
II  describes  the  stationary  model  and  the  approach,  while  Chapter  III  describes  the 
maneuver  model  path.  Once  the  interference  environment  is  characterized,  STAP  can 
be  applied  in  Chapter  IV.  The  stationary  model,  the  foundation  of  the  maneuver 
model  developed  in  Chapter  III,  is  presented  next. 
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III.  Platform  Maneuver  Model  for  Side  Looking  Arrays 


The  popular  model  described  in  Chapter  II  assumed  the  aircraft  platform  was  per¬ 
fectly  stationary  throughout  the  coherent  processing  interval  (CPI).  Abandon¬ 
ing  that  assumption  and  allowing  platform  maneuver  creates  a  more  realistic  model 
and  is  the  focus  of  this  research.  An  important  point  concerning  the  meaning  of 
“maneuver”  in  this  research.  Range  walk  refers  to  the  aircraft  moving  forward  in  the 
direction  of  travel  during  the  CPI.  Range  walk  decorrelation  effects  have  been  de¬ 
rived  and  simulated  in  [5].  The  focus  of  this  research  is  platform  rotation,  not  range 
walk.  Therefore,  the  aircraft  remains  fixed  in  space  while  the  platform  rotates  for  this 
model. 

The  previous  model  of  [3, 10]  serves  as  the  maneuver  model  foundation.  Devel¬ 
opment  of  the  platform  maneuver  STAR  model  begins  with  the  antenna  and  physical 
geometry.  A  mathematical  framework  incorporating  platform  maneuver  about  any 
axis  is  developed.  Once  the  framework  is  constructed,  the  rotation  behavior  of  the 
platform  itself  is  defined.  The  result  is  propagated  through  the  development  of  new 
signal,  target,  and  clutter  models.  The  thermal  noise  model  is  unaffected  by  platform 
maneuver  and  remains  unchanged.  Major  deviations  from  the  previous  model  include: 
a  single  location  now  has  M  different  (6*,0)  pairs  and  M  k  unit  vectors,  the  aircraft 
velocity  vector  changes  each  pulse,  targets  and  clutter  are  illuminated  differently  each 
pulse,  and  target  and  clutter  Doppler  frequencies  change  each  pulse.  Additionally, 
modification  of  the  the  stationary  model  formatting  is  required.  All  steering  vectors 
in  the  maneuver  model  are  MNP  x  1,  versus  Mxl,A^xl,orPxlin  the  stationary 
model.  Maneuver  model  space-time  snapshots  are  created  with  the  Hadamard  prod¬ 
uct  (©),  versus  the  Kronecker  product  (®)  in  the  stationary  model.  The  foundation 
of  these  differences  begins  with  the  geometry  and  rotation  behavior. 

3.1  Geometry  and  Rotation 

Consider  the  side- looking  (SL)  antenna  array  mounted  on  an  aircraft,  as  de¬ 
picted  in  Figures  3.1-  3.2  and  discussed  in  Chapter  II.  Recall  the  unit  vector  from 
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Figure  3.1:  Antenna  array  geometry.  The  antenna  array  is  oriented  in  the  x-z  plane. 
Its  reference  element,  defined  as  the  upper- left  most  element,  serves  as  the  origin  of 
the  cartesian  and  radar  coordinate  systems.  The  elements  are  uniformly  spaced  with 
inner-element  spacings  of  dx  in  the  x  direction  and  dz  in  the  2;  direction.  The  N 
azimuth  channels  are  numbered  0  to  A  —  1 ,  and  the  P  elevation  channels  0  to  P  —  1 . 

the  reference  element  to  the  element  is  [3] 

dnp  =  -nd^-k  -  pdzZ.  (3.1) 


!  I  The  antenna  array  always  defines  the  reference  coordinate  system  for 
this  research.  Any  platform  rotation  equates  to  a  rotation  of  the  reference 
coordinate  system. 
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Reference  Element 


Figure  3.2:  Radar  platform  geometry.  The  array  is  mounted  on  an  aircraft  in  a 

side-looking  (SL)  configuration  with  Foresight  defined  in  the  y  direction.  The  aircraft 
travels  at  speed  Va  in  the  x  direction.  The  vector  k  is  the  unit  vector  from  the  reference 
element  to  the  target. 

The  {x,y,z)  coordinates  of  a  target  or  clutter  patch  at  elevation  angle  9,  azimuth 
angle  0,  and  range  Rt  are 


Rt  cos  9  sin  (p, 

(3.2) 

Rt  cos  9  cos  (p, 

(3.3) 

and 


2;  =  Rt  sin  9. 


(3.4) 
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The  elevation  and  azimuth  angles  can  be  expressed  as 


The  negative  sign  in  Equation  (3.5)  ensures  the  elevation  angle  to  the  ground  is 
negative,  as  expected  in  radar  coordinates.  Figure  3.2  never  changed  in  the  stationary 
model  and  introduction  of  maneuver  at  this  point  serves  as  the  start  of  the  maneuver 
model. 

The  platform  is  allowed  to  rotate  about  the  reference  element  about  any  axis. 
Rotations  about  the  x,y,  and  .2  axes  depicted  in  Figure  3.3  correspond  to  aircraft  roll, 
pitch,  and  yaw  respectively.  Consider  an  arbitrary  point  located  at  {6,  0)  in  Figure 
3.2.  The  antenna  then  rotates  about  any  axis.  That  point  at  (6*,  0)  still  resides  in 
the  same  physical  space,  but  its  position  relative  to  the  antenna  has  changed.  The 
overall  goal  is  accurate  characterization  of  the  clutter  and  target.  The  snapshots 
depend  on  the  scenario  geometry  and  the  resulting  radar  angles  (and  Doppler  values) 
to  the  location.  Therefore,  the  (6*,  (p)  location’s  new  position  relative  to  the  rotated 
antenna  must  be  found.  A  mathematical  framework  is  needed  to  describe  the  radar 
coordinates  to  any  location,  with  any  amount  of  platform  pitch,  roll,  or  yaw.  The 
mathematical  framework  is  constructed  with  rotation  matrices. 

A  general  linear  rotation  about  an  axis  can  be  described  with  a  rotation  matrix. 
Rotations  about  all  three  {x,y,z)  axes  can  be  described  using  Euler’s  rotation  theorem. 
No  universal  convention  exists  defining  rotation  direction,  i.e.  clockwise  (CW)  or 
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Figure  3.3:  Rotation  illustration.  The  antenna  rotates  abont  the  x,  y,  or  z  axis 

corresponding  to  roll  (middle  picture),  pitch  (top  picture),  and  yaw  (bottom  picture) 
respectively.  Rotation  occurs  about  the  reference  element  and  the  antenna  array  re¬ 
sides  in  the  x'-z'  plane.  The  reference  coordinate  system  therefore  rotates,  designated 
by  the  primed  notation.  Rotation  can  occur  in  a  clockwise  or  counter-clockwise  direc¬ 
tion,  althongh  only  clockwise  rotation  is  depicted  above.  The  subscript  corresponds 
to  the  type  of  rotation.  Snbscript  “p”  is  pitch,  “r”  is  roll,  and  “y”  is  yaw. 
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Figure  3.4:  Pitch  up  geometry.  The  array  reference  element  serves  as  the  origin  of 
the  x-y-z  and  x'-y'-z'  cartesian  coordinate  systems.  The  coordinate  system  rotates 
degrees  in  a  CW  direction  about  the  reference  element  and  y-axis,  corresponding  to 
the  aircraft’s  nose  rotating  upward.  The  rotated  coordinates  {x'-^^,y[,z[)  to  the  point 
(xi^yi^Zi)  are  then  found. 

counter-clockwise  (CCW).  The  convention  used  for  this  research  is  based  on  the  pilot’s 
perspective. 


!  I  ‘Pitch  up’  refers  to  a  CW  rotation  about  the  y-ax\s,  i.e.  the  aircraft’s 
nose  rotates  skyward.  ‘Roll  right’  refers  to  a  CCW  rotation  about  the  x-axis, 
i.e.  the  aircraft’s  right  wing  dips  and  left  wing  rises.  ‘Yaw  right’  refers  to  a 
CW  rotation  about  the  2:-axis,  i.e.  the  aircraft’s  nose  rotates  right. 

Consider  the  case  of  pitch  up,  depicted  in  Figure  3.4.  The  original  target  is 
located  at  (xi,zi).  The  platform  rotates  CW  about  the  y-axis  at  the  origin  of  the 
x-y-z  coordinate  system  (the  reference  element).  The  original  target  has  not  moved 
and  exists  at  the  exact  same  location  in  physical  space.  However,  the  relative  frame 
of  reference  defined  by  the  antenna  array’s  location  has  changed.  Inspection  of  Figure 
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3.4  yields  the  new  coordinates  of  the  target  in  the  rotated  x'-y'-z'  coordinate  system, 


x\  =  Xi  cos  ipp  +  Zi  sin  'ipp, 

y'l  =  yi, 

—  —Xi  sin  +  Zi  cos  V'p- 


(3.8) 

(3.9) 
(3.10) 


The  relationship  can  be  expressed  in  matrix  form, 


x'l 

cos  Ipp  0  sin  ipp 

Xi 

y'l 

= 

0  1  0 

yi 

4 

—  sin  ipp  0  cos  ipp 

Zl 

The  3  X  3  rotation  matrix  describes  the  location  of  the  target  with  respect  to  the 
y-axis  rotated  platform  frame  of  reference.  The  matrix  is  dubbed  P“  with  subscript 
p  referring  to  rotation  type  pitch  and  superscript  u  to  rotation  direction  up, 


PW 


cos  ipp  0  sin  ipp 
0  1  0 
—  sin  ipp  0  cos  -ipp 


The  next  case  considered  is  roll  right,  depicted  in  Figure  3.5.  The  platform  rotates 
ipr  degrees  about  the  reference  element  and  x-axis.  Inspection  of  Figure  3.5  reveals 
the  relationship. 


x\ 

=  Xi, 

(3.11) 

y'l 

=  yi  cos  Ipr  +  Z\  sin  Ipr-, 

(3.12) 

A 

=  —y  sin  Ipr  +  z  cos  Ipr- 

(3.13) 
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Figure  3.5:  Roll  right  geometry.  The  array  reference  element  serves  as  the  origin  of 
the  x-y-z  and  x'-y'-z'  cartesian  coordinate  systems.  The  coordinate  systems  rotates 
il^r  degrees  in  a  CCW  direction  about  the  reference  element  and  x-axis,  corresponding 
to  the  aircraft’s  right  wing  dipping  downward.  The  rotated  coordinates  {x'l^y'i^z'i)  to 
the  point  {xi,yi,Zi)  are  then  found. 

The  relationship  in  matrix  form  is, 


1 

o 

o 

Xi 

y'l 

= 

0  cos  Xpr  sin  Ipr 

yi 

. 

0  —  sin  ipr  cos  'ipp 

This  3  X  3  rotation  matrix  describes  the  location  of  the  target  with  respect  to  the 
x-axis  rotated  platform  frame  of  reference.  The  matrix  is  dubbed  with  subscript 
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Figure  3.6:  Yaw  right  geometry.  The  array  reference  element  serves  as  the  origin  of 
the  x-y-z  and  x'-y'-z'  cartesian  coordinate  systems.  The  coordinate  systems  rotates 
ijjy  degrees  in  a  CW  direction  about  the  reference  element  and  .s-axis,  corresponding 
to  the  aircraft’s  nose  rotating  to  the  right.  The  axis  label  “y”  and  rotation  type 
designated  by  subscript  “y”  in  ^|;y  are  nnrelated.  The  notation  ambiguity  is  cleared 
with  the  context  of  nse.  The  rotated  coordinates  {x[,y[,z[)  to  the  point  {xi,yi,zi)  are 
then  found. 

r  referring  to  rotation  type  roll  and  snperscript  r  to  rotation  direction  right, 


PMv) 


1  0  0 

0  cos  'Ipr  sin 
0  —  sin  cos 


Rotation  about  the  y  (pitch)  and  x  (roll)  axes  have  been  defined,  leaving  rotation 
abont  the  z-axis  (yaw).  Figure  3.6  depicts  the  platform  yawing  right  a  total  of  ipy 
degrees.  Inspection  of  Fignre  3.6  reveals  the  relationship, 


x\  =  Xi  cos  xpy  —  yi  sin  'ipy, 
y[  =  xi  sin  i)y  +  yi  cos 


(3.14) 

(3.15) 

(3.16) 


Zi  =  Zi- 


And  the  relationship  in  matrix  form  is, 


x[  cos  -ipy  —  sin  'ipy  0  xi 

y[  =  sin 'ipy  cosily  0  yi 
z[  0  0  1  Zi 


This  3  X  3  rotation  matrix  describes  the  location  of  the  target  with  respect  to  the 
2;-axis  rotated  platform  frame  of  reference.  The  matrix  is  dubbed  with  subscript 
r  referring  to  rotation  type  yaw  and  superscript  r  to  rotation  direction  right, 


cos  iy  —  sin  iy  0 
=  siu'^j^  cos  iy  0 

0  0  1 


Realistically,  the  platform  is  not  restricted  to  upward  pitch,  right  roll,  and  right  yaw. 
The  platform  may  maneuver  opposite  each  of  these  respective  directions,  necessitat¬ 
ing  examination  of  downward  pitch,  left  roll,  and  left  yaw.  These  three  maneuver 
scenarios  are  depicted  in  Figures  3. 7-3. 9.  And  the  rotation  matrices  corresponding 
to  pitch  down,  roll  left,  and  yaw  left  are 


cos  ip  0  —  sin  ip 


0  1  0 


sin  ip  0  cos  ip 


1  0  0 

Pr(l/'r)  =  0  cosir  —  Smir 

0  sin  ir  COS  ir 
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Figure  3.7:  Pitch  down  geometry.  The  array  reference  element  serves  as  the  origin 
of  the  x-y-z  and  x'-y'-z'  cartesian  coordinate  systems.  The  coordinate  systems  rotates 
ipp  degrees  in  a  CCW  direction  about  the  reference  element  and  zy-sods,  corresponding 
to  the  aircraft’s  nose  rotating  to  the  right.  The  rotated  coordinates  {x[,y[,z[)  to  the 
point  {xi,yi,zi)  are  then  found. 


and 


?;,(</>,)  = 


cos  Ipy 

—  sin  Ipy 
0 


sin  xpy  0 
cos  xpy  0 
0  1 
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Figure  3.8:  Roll  left  geometry.  The  array  reference  element  serves  as  the  origin  of 
the  x-y-z  and  x'-y'-z'  cartesian  coordinate  systems.  The  coordinate  systems  rotates 
ipr  degrees  in  a  CCW  direction  about  the  reference  element  and  x-axis,  corresponding 
to  the  aircraft’s  left  wing  dipping  and  right  wing  rising.  The  rotated  coordinates 
to  the  point  {xi,yi,zi)  are  then  found. 

respectively.  The  six  rotation  matrices,  P",  P^,  Pj;,  Pj,,  P^,  P^,  enable  description 
of  any  rotation  (CW  or  CCW)  about  any  axis.  Each  rednces  to  an  identity  matrix 
when  ipp  =  tpr  =  P’y  =  0,  as  expected.  Examination  of  pitch  matrices  P“  and  P^ 
reveal  the  matrices  are  nearly  identical  with  the  exception  of  the  sign  of  sin  pjp.  The 
same  similarity  holds  for  roll  matrices  PJ^  and  P^,  and  yaw  matrices  P^  and  P^.  The 
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Figure  3.9:  Yaw  left  geometry.  The  array  reference  element  serves  as  the  origin  of 
the  x-y-z  and  x'-y'-z'  cartesian  coordinate  systems.  The  coordinate  systems  rotates 
ijjy  degrees  in  a  CCW  direction  about  the  reference  element  and  2:-axis,  corresponding 
to  the  aircraft’s  nose  rotating  to  the  left.  The  rotated  coordinates  {x']^,y[,z[)  to  the 
point  {xi,yi,zi)  are  then  found. 

oddness  property  of  the  sin  function  allows  reduction  from  the  six  rotation  matrices 
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to  three:  Pp,  P^,  and  P^, 


and 


COS  A  0 

sin  A 

Pp('0p)  — 

0  1 

0 

— 

sin  A  0 

cos  A 

1 

0 

0 

PriA)  = 

0 

cos  A 

sin  A 

0 

—  sin  ipr 

cosV'r 

cos  ipy  —  sin  xpy  0 
sin  xjjy  cos  xjjy  0 

0  0  1 


!  I  Convention  for  this  research: 

A  positive  ipp  refers  to  a  CW  rotation  about  the  y-axis,  i.e.  ‘Pitch  Up’. 

A  positive  pjr  refers  to  a  CCW  rotation  about  the  a:-axis,  i.e.  ‘Roll  Right’. 
A  positive  pjy  refers  to  CW  rotation  about  the  2:-axis,  i.e.  ‘Yaw  Right’. 

The  complete  rotation  matrix  describing  any  rotation  three  axes  is  found  by 


However,  the  order  of  operations  will  affect  the  rotation  being  modeled.  The  rota¬ 
tion  Pp{'lpp)  X  [Pr{'lpr)  X  ^yii^y)]  is  HOt  the  Same  as  [Pp{'lpp)  X  Pri^r)]  X  Pyi^pr)  or 
\Py{'ipp)  X  Pri^ipy)]  X  Pp{'ipr)-  In  other  words,  the  order  of  successive  rotations  impacts 
the  overall  rotation  and  the  resulting  geometry.  A  simple  example  depicted  in  Fig¬ 
ure  3.10  illustrates  the  concept.  Three  rotations  are  performed  in  different  orders, 
resulting  in  two  different  results.  The  different  resulting  geometries  is  a  consequence 
of  the  approach  taken  in  this  research.  Specifically,  the  rotation  matrices  describe  the 
location  of  point  {x,y,z)  in  the  {x',y',z')  rotated  reference  coordinate  system.  The  as- 
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90  CW  ~X  90  CCW  ~  X 
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90cw~  y 


Figure  3.10:  Rotation  order  impact.  The  axes  undergo  three  rotations:  a  90°  CW 
rotation  about  the  a:- axis,  a  90°  CW  rotation  about  the  y-axis,  and  a  90°  CCW 
rotation  about  the  x-axis.  The  rotations  are  applied  in  different  orders,  resulting  in 
different  ending  geometries. 

sumption  in  this  research  is  the  rotated  reference  coordinate  system  does  not  remain 
fixed.  In  the  maneuver  model,  six  (6)  unique  rotation  orders  exist,  each  with  their  own 
unique  rotation  matrix:  pitch-roll-yaw,  pitch-yaw-roll,  roll-pitch-yaw,  roll-yaw-pitch, 
yaw-roll-pitch,  and  yaw-pitch-roll.  It  is  also  feasible  to  change  rotation  order  during 
the  coherent  processing  interval  (CPI),  i.e.  perform  a  pitch-roll-yaw  rotation  on  one 
pulse,  then  a  yaw- pitch- roll  rotation  on  the  next.  The  focus  of  this  research  is  char¬ 
acterization  of  individual  platform  maneuver  effects.  While  the  rotation  framework 
constructed  allows  any  combination  of  maneuver,  each  will  be  investigated  separately. 

A  brief  summary  is  in  order.  Rotation  matrices  for  aircraft  pitch,  roll,  and  yaw 
were  constructed  using  a  relative  reference  frame  (the  platform).  These  rotation  ma¬ 
trices  produce  the  {xyy[,z[)  coordinates  of  the  original  [xi^yi^Z])  location.  The  overall 
goal  remains  characterization  of  the  target  and  clutter  environment.  The  (a:(,|/(,;2() 
coordinates  are  therefore  transformed  into  their  corresponding  radar  coordinates  (d^, 
0'^).  Analogously  to  Equations  (3.5)  and  (3.6),  the  elevation  angle  9[  and  azimuth 
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angle  (/)[  to  a  target  with  coordinates  at  range  Rt  are  given  by 


sin 


(3.18) 


and 


cos 


-1 


(3.19) 


However,  Equation  (3.19)’s  general  definition  of  azimuth  angle  will  not  suffice  in  the 
platform  maneuver  model.  Equation  (3.19)  is  based  on  the  geometry  of  Eigure  3.2. 
It  provides  the  angle  from  the  axis  on  which  y[  exists,  as  depicted  in  Eigure  3.11. 
This  subtlety  was  a  non-issue  in  the  stationary  model,  since  the  azimuth  angle  to  a 
target  or  clutter  patch  was  constant  and  never  required  recalculation.  A  target  was 
simply  placed  at  an  elevation  and  azimuth  angle,  while  the  ground  plane  was  divided 
into  i  X  k  patches  and  the  azimuth  angles  to  the  patches  were  Nc  linearly  spaced 
angles  from  0-27r.  The  maneuver  model  azimuth  angle  is  therefore  first  calculated 
using  Equation  (3.19),  then  modified  appropriately  depending  on  the  quadrant  the 
original  {xi,yi,zi)  location  resides  in  the  rotated  {x'-^^,y[,z[)  space: 


/ 


</>i  A,  y[>0 

>  0  and  y[  <  Q 

0i-7r  x\,y[<{) 

—(p'l  x\  <  0  and  >  0 


Expressions  for  azimuth  and  elevation  angle  for  pitch,  roll,  can  now  be  derived. 
Recall  {x'l^y'i^z'i)  are  functions  of  rotation  type,  rotation  amount,  and  the  original 
(xi,r/i,.2i)  coordinates.  In  the  case  of  pitch.  Equations  (3.8  -  3.10)  are  substituted  for 
(xi,yi,2;i)  in  Equations  (3.18  -  3.19)  producing. 


(3.20) 
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Figure  3.11:  Azimuth  angle  calculation.  The  azimuth  angle  to  the  loca¬ 

tion  depends  on  the  quadrant  it  resides  in  the  rotated  coordinate  system.  Eqnation 
(3.19)  calcnlates  the  angle  from  the  r/'-axis  in  the  specific  qnadrant.  The  resnlt  mnst 
therefore  be  translated  to  the  corresponding  angle  from  boresight. 

and 


0;  = 


COS 


-1 


Rt  cos  Op 


(3.21) 


Snbstitnting  Eqnations  (3.2  -  3.4)  into  (3.20  -  3.21)  yields  the  elevation  and  azimnth 
angles  to  the  original  {xi,yi,zi)  given  any  amount  of  pitch 


6p  =  sin  ^  (— cos6*sin0sin'^p -|- sin6*cos'^p) , 


(3.22) 
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and 


-1 


cos  9  cos  (j) 


(3.23) 


cos 


Pitch  effect  on  azimuth  and  elevation  angle  is  illustrated  in  Figure  3.12.  An  (x,  y, 
z)  ground  location  was  chosen  in  each  quadrant.  The  degree  of  pitch  was  varied 
from  0^-360'^  corresponding  to  the  platform  making  one  full  rotation  about  its  x-axis. 
While  unrealistic  for  the  airborne  radar  platform,  the  full  rotation  allows  insight  into 
pitch  effects  on  the  angles  to  various  locations.  The  results  match  physical  intuition. 
Consider  the  ground  location  at  boresight  (0  =  0°),  labeled  with  a  diamond  (o)  in  the 
plot.  The  location  is  originally  located  at  —3°  elevation.  Recall  the  definition  of  the 
elevation  angle  of  Equation  (3.18),  stating  the  elevation  angle  is  dependant  on  z' .  As 
the  platform  begins  to  pitch  up  (CW),  the  z'-axis  rotates  CW,  the  z'  coordinate  of  the 
ground  location  becomes  less  negative,  and  Equation  (3.18)  approaches  zero.  When 
the  pitch  rotation  is  90°,  the  ground  location  is  at  z'  =  0  and  the  elevation  angle  to 
the  original  location  is  now  0°.  Continuing  the  rotation  past  90°  places  the  platform 
upsidedown,  and  the  elevation  angle  to  the  original  location  becomes  positive.  At 
180°  pitch,  the  platform  is  completely  inverted  and  the  elevation  angle  to  the  location 
is  +3°.  Similar  basic  exercises  for  the  other  locations  plotted  and  azimuth  angles  show 
the  results  match  physical  intuition. 

Next,  roll  effects  on  elevation  and  azimuth  angles  are  found  by  substituting 
Equations  (3.12  -  3.13)  for  (xi,yi,Zi)  in  Equations  (3.18  -  3.19)  producing. 


6'j.  =  sin  ^  (— cos0cos(/)sin'i/ir  +  sin^cos'^r) 


(3.24) 


and 


(3.25) 


Roll  effects  on  azimuth  and  elevation  angle  are  illustrated  in  Figure  3.13.  Once  again 
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Figure  3.12:  Pitch  effect  on  azimuth  and  elevation  angles.  Five  different  azimnth 

locations  on  a  single  range  ring  are  plotted.  The  top  figure  shows  the  pitch  effect  on 
elevation  angle  to  various  ground  locations  at  a  range  of  60  km.  The  bottom  figure 
shows  the  effect  on  azimuth  angle  to  those  same  locations.  The  sharp  jumps  to/from 
±180°  correspond  to  the  location  passing  throngh  radar  the  Foresight  {y')  axis. 
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Figure  3.13:  Roll  effect  on  azimuth  and  elevation  angles.  Five  different  azimuth 

locations  on  a  single  range  ring  are  plotted.  The  top  figure  shows  the  roll  effect  on 
elevation  angle  to  various  ground  locations  at  a  range  of  60  km.  The  bottom  figure 
shows  the  effect  on  azimuth  angle  to  those  same  locations.  The  sharp  jumps  to/from 
±90°  correspond  to  the  location  passing  to/from  the  front/back  side  on  the  array. 
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physical  intuition  and  the  results  match.  Consider  the  ground  location  at  boresight 
(0  =  0),  labeled  with  a  diamond  (o)  in  the  plot.  Elevation  angle  changes  linearly  with 
'ipr  for  the  location  0  =  0,  and  azimuth  angle  to  0  =  0  is  constant  despite  the  roll. 


And  hnally,  the  elevation  and  azimuth  angles  for  yaw  are  found  by  substituting 
Equations  (3.15  -  3.16)  for  (xi,r/i,.2i)  in  Equations  (3.18  -  3.19)  producing. 


6'y  =  sin  ^  (sind)  =  d, 


(3.26) 


and 


(3.27) 


The  result  of  Equation  (3.26)  states  the  elevation  angle  to  any  location  is  unchanged, 
despite  any  amount  of  platform  yaw.  The  result  makes  physical  sense  since  yaw 
corresponds  to  a  rotation  of  the  x-y  plane.  The  elevation  angle  9  is  measured  from 
the  x-y  plane  to  the  target,  as  depicted  in  Figure  3.2.  The  rotated  x-y  plane  remains  in 
the  original  plane  for  any  values  of  yaw  and  more  importantly  the  point  from  which 
9  is  measured  does  not  move.  Therefore,  the  elevation  angle  remains  unchanged. 
Substituting  9y  =  9  into  the  expression  for  0(^  (Equation  3.27)  gives 


y 


cos 


(3.28) 


The  cos  9  terms  cancel,  and  the  compound  angle  identity 


sin  A  sin  B  +  cos  A  cos  B  =  cos(A  —  B) 


(3.29) 


produces 


y 


cos  ^  (cos(0  —  ipy)) 


(3.30) 

(3.31) 


0-0, 
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Equation  (3.31)  demonstrates  the  azimuth  angle  in  the  rotated  coordinate  system, 
is  a  linear  function  of  the  yaw  amount  i/^y.  Yaw  effects  are  illustrated  in  Figure 
3.14.  The  unchanging  elevation  angle  and  linearity  of  the  azimuth  angle  are  shown. 

In  this  section,  a  mathematical  framework  enabling  description  of  the  radar  an¬ 
gles  to  any  location  with  any  amount  of  pitch,  roll,  or  yaw  was  constructed.  With 
this  mathematical  framework  in  place,  the  logical  progression  is  characterization  of  the 
platform  rotation  itself.  Equations  (3.22  -  3.27)  were  defined  for  a  total  degree  of  rota¬ 
tion,  but  the  platform  maneuvers  continuously  throughout  the  CPI.  That  movement 
must  be  discretized  and  the  total  degree  of  rotation  at  each  discretization  quantified 
in  order  to  apply  the  framework  developed  in  this  section. 


3.2  Rotation  Behavior 

The  first  step  in  describing  the  rotation  behavior  is  defining  the  time  space. 
The  radar  transmits  an  uncompressed  pulse  of  width  r  seconds  then  samples  the 
environment  at  the  center  of  each  range  cell,  as  depicted  in  Figure  3.15.  The  radar 
begins  transmitting  at  time  t  =  0,  ends  transmitting  at  time  t  —  r,  then  samples  L 
range  cells  at  r  intervals. 


!  For  this  development,  the  platform  rotates  continuously.  The  platform 
orientation  on  transmit  is  therefore  different  than  the  orientation  on  receive. 


The  process  is  repeated  M  times,  once  for  each  pulse  in  the  CPI.  The  entire 
time  span  of  the  CPI  is  therefore  M{L  +  1)  samples,  beginning  at  the  transmit  time 
for  the  zero*^  pulse  (|)  and  ending  at  the  sampling  time  for  the  range  cell  on  the 
(M  —  1)*^  pulse. 


t  =  T 


1 

2 


:  1  :  (M(L  +  1)) 


T 


(3.32) 


The  transmit  and  receive  times  are  subsets  of  the  entire  time  span.  The  first  time 
sample  |  corresponds  to  the  first  transmitted  pulse.  The  next  pulse  is  transmitted 
L  -|-  1  samples  later,  after  the  radar  has  sampled  the  L  range  cells.  The  third  pulse  is 
transmitted  L  -|- 1  samples  after  the  second,  or  2K  -|-  2  after  the  first.  The  M  transmit 
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Figure  3.14:  Yaw  effect  on  azimuth  and  elevation  angles.  Five  different  azimuth 

locations  on  a  single  range  ring  are  plotted.  The  top  ffgure  shows  the  yaw  effect 
on  elevation  angle  to  various  ground  locations  at  a  range  of  60  km.  As  shown  in 
Equation  3.26,  yaw  has  no  effect  on  elevation  angle.  The  bottom  figure  shows  the 
effect  on  azimuth  angle  to  those  same  locations.  As  shown  in  Equation  3.31,  the 
azimuth  angle  changes  linearly  with  the  amount  of  yaw.  The  sharp  jumps  at  various 
locations  correspond  to  each  location  passing  through  the  radar  Foresight  axis. 
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Figure  3.15:  Maneuver  model  pulse  train  and  time  span.  The  CPI  begins  at  t  =  0 
with  the  platform  in  its  original  position.  Platform  rotation  begins  immediately  at 
t  =  0.  The  platform  rotation  is  discretized  into  r  intervals.  Transmit  times,  noted 
with  diamonds  on  the  pnlses,  are  designated  txmit  and  the  corresponding  azimuth  and 
elevation  angle  to  the  location  (d,0)  are  designated  Receive  times,  noted 

with  diamonds  on  range  cells,  are  designated  trcv  and  the  corresponding  azimuth  and 
elevation  angle  to  the  location  (d,^)  are  designated  There  are  M  values  for 

each  of  these  quantities,  corresponding  to  the  M  pulses  in  the  CPI.  Range  cells,  k, 
are  nnmbered  1  to  L. 

times  are  therefore  the  subset  of  the  entire  time  span  given  by, 


fymit  —  t(l  :  T  +  1  :  M(L  +  1)) 


(3.33) 


The  receive  times  depend  on  the  range  cell  under  test  (RUT).  A  CPI  contains  M 
different  returns  sampled  at  M  different  times  corresponding  to  the  same  range  Rt. 
The  total  number  of  range  cells  is  defined  as 


L  = 


(3.34) 


where  [  J  represents  the  flooring  function.  The  range  bin  of  interest,  k,  corresponding 
to  the  time  sample  at  range  Rt  is 


K  = 


RtL 

Ru  ’ 


(3.35) 
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where  is  the  radar’s  unambiguous  range, 


(3.36) 


The  first  receive  time  corresponds  to  the  k  +  1  sample,  as  depicted  in  Figure  3.15. 
The  next  receive  time  is  L  +  1  samples  later.  And  the  third  receive  time  is  L  +  1 
samples  after  the  second,  or  2L  +  2  after  the  first.  The  M  receive  times  are  therefore 
a  subset  of  the  entire  time  span  given  by 


trcv  =  t(K  +  1  :  L  +  1  :  M(L  +  1)). 


(3.37) 


Given  complete  characterization  of  the  CPI  time  span,  the  platform  orientation  at 
each  transmit  and  receive  time  must  be  defined.  Pitch,  roll,  and  yaw  rates  (degrees 
per  second)  are  defined  as  Pp,  pr,  and  Py  respectively.  These  rotations  rates  are 
assumed  constant  throughout  the  CPI. 


!  I  As  previously  mentioned,  the  focus  of  this  research  is  the  effect  of 
individual  rotations.  While  the  development  of  each  rotation  to  follow  will 
mirror  one  another,  only  one  rotation  type  will  be  implemented  at  a  time. 
Equations  (3.38  -  3.40)  yield  M  linear  rotation  angles  based  on  the  constant 
rotation  rates.  If  modeling  combinations  of  maneuver  effects  were  of  interest, 
a  deviation  from  the  constant  rotation  rate  could  be  employed,  i.e.  platform 
rolls  right  for  several  time  samples,  then  yaws  up,  then  rolls  left.  The  constant 
rotation  rate  limits  the  modeling  capability  to  linear  combinations  (pitch-roll- 
yaw,  pitch-yaw-roll,  roll-pitch-yaw,  roll-yaw-pitch,  yaw- pitch- roll,  or  yaw-roll- 
pitch)  of  a  set  number  of  degrees  each  time  sample. 

At  any  given  time  the  total  degree  of  rotation  is  the  product  of  rotation  rate 
and  time.  The  total  degree  of  rotation  at  each  instance  of  transmit  are  therefore  the 
rotation  rates  times  the  transmit  times. 


(3.38) 

(3.39) 
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and 


'4^y  —  Pj/txmit-  (3.40) 

And  the  total  degree  of  rotation  at  each  instance  of  receive  are  the  rotation  rates  time 
the  receive  times, 


=  pptrcv{Rt),  (3.41) 

■07  =  PrircviRt),  (3.42) 


and 


0y  Py^rcv(^Rt}  ■ 


(3.43) 


!  I  All  angles  corresponding  to  array  orientation  upon  transmit  are  desig¬ 
nated  with  superscript  tx.  All  angles  corresponding  to  array  orientation  upon 
receive  are  designated  with  superscript  rv.  The  subscript  indicates  the  type 
of  rotation. 

The  M  elevation  angles  to  the  location  (d,0)  on  transmit  are  therefore  the 
rotation  amounts  (Equations  (3.38-3.40))  evaluated  by  the  expressions  for  elevation 
angle  (Equations  (3.22-3.26)): 


=  sin  ^  (— cos6*sin0sin0p^ -h  sin6*cos0p^)  ,  (3.44) 


07  =  sin  ^  (—cos d cos 0 sin 07  +  sill ^ cos 07)  )  (3.45) 


or 


6>7  =  eiM,  (3.46) 

depending  on  the  rotation  being  modeled.  The  1m  represents  a  M  x  1  vector  of  ones, 
and  therefore  6*OlMisaMxl  vector  with  M  values  of  the  original  6  (yaw  does  not 
affect  elevation  angle). 
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!  A  element  by  element  division  operator  must  be  introduced  to  proceed. 
For  this  research,  element  by  element  division  is  represented  with  the  □  in 
the  denominator  of  an  expression.  The  □  precedes  the  vector  quantity  that 
performs  the  element  by  element  division. 

The  M  azimuth  angles  to  location  {9,(/))  on  transmit  are: 


(pf  =  cos 


-1 


cos  6  cos  (p 
□  cos0‘^ 


(3.47) 


€  = 


cos 


-1 


cos  9  cos  (p  cos  -0*^  +  sin  9  sin  ip. 


'  tx 


□  cos  6 


tx 


or 


€  = 


cos 


-1 


cos  9  sin  (/)  cos 


tx 


cos  9  cos  p  cos  ipy 


tx 


□  cos  0 


tx 


(3.48) 


(3.49) 


depending  on  the  rotation  being  modeled.  And  the  M  elevation  and  azimuth  angles 
on  receive  are: 


0)7  =  sin  ^  (— cos 0 sine/) sin 0)7  +  sill ^ cos 0)7)  )  (3.50) 


0"  =  sin  ^  (— cos0cos(/)sin0)'' +  sin0cos0)'') ,  (3.51) 


or 


0 


rv 

y 


91m, 


and 


0 


rv 

P 


COS 


cos  9  cos  p\ 
□  cos0"  )  ’ 


(3.52) 


(3.53) 
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4>r 


COS 


(3.54) 


-1 


COS  6  cos  (p  cos  'tj/pp  +  sin  6  sin  ip 


□  cos  91 


or 


,  /  cos  9  sin  (h  cos  iblj  +  cos  9  cos  6  cos  iblj  \ 

<f7  -  --  ( -  □cos.-  ’)  ’ 

depending  on  the  rotation  being  modeled. 

The  mathematical  framework  characterizing  aircraft  pitch,  roll,  or  yaw  in  either 
CW  or  CCW  directions  was  established  in  Section  3.1.  The  rotation  behavior  was 
characterized  in  this  section.  The  net  result  is  locations  (targets  and  clutter  patches) 
previously  dehned  by  a  single  {9,(f))  pair  now  are  defined  by  M  (0*^,  (fp^)  pairs  on 
transmit  and  M  (0-,  0-)  pairs  on  receive.  The  objective  remains  accurate  charac¬ 
terization  of  the  target  and  the  clutter  environment.  The  stationary  model  framework 
discussed  in  Chapter  II  serves  as  the  foundation  for  that  characterization,  but  requires 
modification.  Components  of  the  stationary  model  requiring  modification  include  the 
signal  model,  formatting,  target  model,  and  clutter  model. 

3.3  Platform  Maneuver  Signal  Model 

While  the  orientation  of  the  platform  on  transmit  changes  on  a  pulse  to  pulse 
basis,  the  transmit  signal  itself  is  unchanged  and  remains  [3], 

s{t)  =  (3.56) 

However,  modification  of  the  receive  signal  model  is  required.  Quantities  that  were 
time  independent  in  the  stationary  model  now  have  time  dependencies.  Derivation  of 
the  platform  maneuver  receive  signal  model  begins  with  the  original  stationary  model 
receive  waveform  defined  by  [3] 

Snp{t)  =  a^u(t  -  (3.57) 
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where  the  time  delay  Tnp  is  the  sum  of  the  round  trip  time  from  the  array  face  to  the 
target,  Tr,  and  the  differential  time  from  the  reference  element  to  each  individual 
npth  element  [3], 

Tnp  =  Tr  +  r'p.  (3.58) 

The  narrowband  assumption, 


T  > 


c 


(3.59) 


remains  and  is  unaffected  by  platform  maneuver.  The  physical  dimensions  of  the 
array  and  time  duration  of  the  signal  have  not  changed.  Each  element  receives  the 
signal  at  the  same  time  under  the  assumption,  and  the  signal  is  rewritten  as 


Snp{t)  =  a^u(t  -  (3.60) 


Platform  maneuver  equates  to  time-varying  array  orientation  and  velocity  vector. 
Therefore,  the  differential  time  delay  t!^p  and  Doppler  frequency  fd  now  have  time 
dependencies.  Additionally,  the  receive  amplitude  is  no  longer  assumed  constant  for 
the  CPI.  The  environment  is  illuminated  differently  as  the  platform  rotates  and  the 
receive  amplitude  therefore  requires  a  time  dependence  as  well.  The  total  time  delay 
becomes 

fnpit)  =TR  +  f^pY).  (3.61) 

And  the  differential  time  delay  with  time  dependance  is  given  by, 


np 


C  ’ 

—ndx  cos  sin  —  pd^  sin 

c 


(3.62) 

(3.63) 


where  the  radar  angles  and  platform  orientation  are  both  now  time  dependant.  The 
time-variant  spatial  frequencies,  describing  the  phase  behavior  of  the  return  signal 
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across  element,  are  given  by 


and 


^xit)  = 


dx  cos  sin 


^z{t)  = 


dz  sin 

a:  ' 


(3.64) 


(3.65) 


Incorporating  these  time  dependencies  into  the  signal  model  yields 

(3.66) 

The  phase  terms  is  not  a  function  of  time  and  therefore  rolled  into  the 

random  phase  term  e^‘^, 

Snp{t)  =  arit)u{t  -  (3.57) 


Referring  to  Equation  3.63,  fdit)^npi^)  ~  0  because  the  element  spacing  is  on  the  order 
of  several  meters  while  the  propagation  speed  is  the  speed  of  light  (on  the  order  of  10® 
meters).  Expected  Doppler  frequencies  are  not  large  enough  to  make  the  fa^npit)  ~  0 
approximation  invalid.  This  reduction  is  not  performed  on  the  foT^pit)  term  because 
the  transmit  frequency  /„  can  be  on  the  same  order  of  magnitude  as  c.  Approximating 
e^fd(t)Tnp(t)  as  e°  =  1  leaves  the  return  signal  given  by 

Snpit)  =  ar{t)u{t  -  r«)e^'2’^(A+A(0)ie-i2^Ah)Tflgi2^M4t)+p,?.(t))gjV_ 


Each  element,  with  its  own  receiver  chain  depicted  in  Eigure  2.3,  strips  the  carrier 
frequency  fo  to  an  intermediate  frequency  assumed  baseband: 

Snp{t)  =  ar{t)u{t  -  r«)e^'2’^A(0te-j’2’^AWTRgi27r(n,?,(t)+p,?,(t))gjV_ 
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Platform  maneuver  has  no  effect  on  the  Doppler  tolerance  condition  of  [5] 


n,.  <  .  (3.70) 

The  time- frequency  autocorrelation  function  is  therefore  assumed  unity  and  the  match 
filtered  signal  is  given  by 


(3.71) 


The  radar  samples  the  environment  at  the  RUT  M  times  during  the  CPI.  The  round 
trip  delay  time  for  the  pulse  is  given  by  Tr  -|-  mTj..  The  match  filtered  signal  at 
pulse  m  is  therefore 


X 


mnp 


M-1 

Q_^(.^)gjVgi27r(ni?a:(m)+pi?2(m))  ^j2'n-fa{i)(Ti{+iTr) 

i=0 


(3.72) 


The  Tr  term  can  be  extracted  from  the  summation  and  rolled  into  the  random  phase 
And  the  term  fa  x  is  the  definition  of  normalized  Doppler  frequency  O.  The 
final  match  filtered  signal  for  pulse  m  is  given  by 


m 

_  ^^^^^gi27r(ni?a:(»ra)+pi?z(m))  gi27ra>(i)^ 

2=0 


(3.73) 


where  and  were  rolled  into  the  complex  amplitude  at(m). 

Reduction  to  the  stationary  model  waveform  is  achieved  when  t?a;(m)  = 

'&z{'n^)  =  ^zi  ^(0  =  Citim)  =  CKf. 

Several  differences  are  apparent  when  comparing  the  maneuver  and  stationary 
signal  models.  The  spatial  frequencies  now  depend  on  the  platform  orientation  at 
pulse  m.  The  phase  due  to  Doppler  is  no  longer  a  linear  function  of  a  single  a),  as  the 
Doppler  from  each  pulse  depends  on  platform  orientation  at  that  pulse.  The  return 
amplitude  is  no  longer  constant  throughout  the  CPI.  The  environment  is  illuminated 
differently  each  pulse  as  the  platform  maneuvers.  In  order  to  accurately  characterize 
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the  target  and  clutter,  space-time  snapshots  must  be  created.  However,  these  signal 
model  differences  created  by  platform  maneuver  necessitate  a  modification  of  the  sta¬ 
tionary  model  framework.  The  time  dependencies  of  spatial  and  Doppler  frequencies, 
as  well  as  amplitude,  does  not  “fit”  into  the  previous  framework.  The  next  section 
creates  a  new  framework  to  create  accurate  space-time  snapshots  and  model  the  target 
and  interference  environment. 

3.4  Formatting 

The  discussion  to  follow  describes  maneuver  steering  vector  and  maneuver  space- 
time  snapshot  construction  to  a  single  location  originally  at  (0,0)  in  the  stationary 
platform  scenario.  The  format  developed  here  is  then  applied  to  the  modified  target 
(Section  3.5)  and  clutter  (Section  3.6)  models.  A  target  or  clutter  patch  was  previously 
located  at  one  9  and  one  0  during  a  single  CPI.  That  single  (0,  0)  translated  to  single 
valued  spatial  frequencies,  k  unit  vector  to  (0,  0),  relative  velocity,  Doppler  frequency, 
platform  orientation,  antenna  pattern  to  (0,  0),  and  return  amplitude.  The  spatial 
steering  vectors  were  N  x  1  and  P  x  1  column  vectors,  while  the  temporal  steering 
vector  was  M  x  1.  The  MNP  x  1  steering  vector  to  (0,  0)  was  created  with  the 
Kronecker  product  (®), 


Vf(0,  0)  =  e(d^)  ®  b((x;)  ®  a(d3,). 

And  the  complete  space-time  snapshot  was 


(3.74) 


Xt(0,0)  =  tttvt.  (3.75) 

Incorporating  platform  maneuver  creates  M  different  values  of  the  aforementioned 
quantities  for  each  (0,0)  location,  necessitating  modification  of  the  Chapter  II  steer¬ 
ing  vector  framework.  Capturing  the  MNP  x  1  space-time  snapshot  describing  the 
per  pulse  per  element  signal  return  remains  the  goal,  but  a  new  steering  vector  and 
amplitude  framework  is  required.  The  steering  vectors  will  no  longer  be  separable  as 
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Figure  3.16:  Maneuver  model  formatting  block  diagram.  The  following  depiction 

applies  to  a  single  location  located  at  (0,0)  in  the  original  coordinate  system.  The 
location  is  transformed  into  M  transmit  and  M  receive  coordinates  using  the  rotation 
matrices  developed  in  Section  3.1.  The  M  pairs  are  then  transformed  into  radar 
coordinates.  The  M  transmit  and  M  receive  angles  to  location  (0,0)  translate  to 
M  different  array  gains,  element  gains,  element  location  nnit  vectors,  nnit  vectors 
to  (0,0),  and  aircraft  velocity  vectors.  M  x  1  amplitudes,  spatial  steering  vectors, 
and  temporal  steering  vectors  result.  These  are  then  recast  into  MNP  x  1  spatial- 
temporal  steering  vectors.  The  snapshot  is  given  by  the  Hadamard  product  of  the 
four  MNP  X  1  vectors. 

spatial  and  temporal  in  the  strictest  sense,  because  each  will  contain  spatial  and  tem¬ 
poral  phase  information.  Figure  3.16  offers  a  visualization  of  the  formatting  approach 
for  the  platform  maneuver  model. 

Section  3.1  derived  the  rotation  matrices  for  pitch,  yaw,  and  roll.  Section  3.2  de¬ 
scribed  the  discretized  platform  maneuver,  resulting  in  differing  platform  orientations 
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on  transmit  and  receive  for  each  pulse.  The  first  maneuver  model  formatting  step  is 
applying  two  (one  for  transmit  and  one  for  receive)  3  x  3  x  M  rotation  matrices, 
and  P",  to  the  original  {x,y,z)  coordinates  of  the  {0,(f))  location  as  depicted  in  Figure 
3.16.  P*^  and  P'^''  are  filled  with  rotation  values  based  on  the  constant  rate  of  rota¬ 
tion  throughout  the  CPI  as  defined  in  Section  3.2.  Equations  (3.18-3.19)  transform 
the  resulting  M  ,z^^)  and  {x^^' ,y^'' ,z’^'')  values  into  M  radar  coordinate  values 

(d‘^,(/)*^)  and  (d",(/)’’'').  The  M  radar  coordinate  values  to  (0,0)  are  cast  into  M  x  1 
column  vectors  0*^, <0*^,0",  and  0". 

Recall  the  k  unit  vector  describes  the  direction  to  the  (0,  0)  location  referenced 
from  the  array  orientation  in  the  rotated  frame.  A  different  target  unit  vector  therefore 
exists  for  each  platform  orientation.  The  platform  orientations  on  receive  provide  the 
needed  phase  information  and  Doppler  information  of  the  return  signal.  The  M  x  1 
column  vector  of  unit  vectors  to  (0,  0)  on  receive  are  designated  and  given  by 

k’''"  =  [cos  0"  ©  sin  -|-  [cos  0^'"  0  cos  +  [sin  (3.76) 

The  array  reference  element  defines  the  reference  frame  origin  at  any  orientation,  and 
the  element  locations  therefore  remain  in  the  x-z  reference  plane.  There  are  now  M 
x-z  reference  planes,  corresponding  to  M  platform  orientations  on  receive.  The  M 
unit  vectors  from  the  reference  element  to  the  element  are  thus  given  by 

^np  =  -  pdzzh"'.  (3.77) 

,  ,  ^  ^  rv 

The  newly  defined  6  ,(j)  0",  and  and  resulting  dimensionality  drive 

the  creation  of  MNP  x  1  spatial-temporal  steering  vectors. 

3.4.1  Spatial- Temporal  Steering  Vectors  with  Motion.  The  differential  time 
delay,  rfp,  describes  the  return  signal  time  difference  of  arrival  between  elements. 
That  time  difference  of  arrival  depends  on  the  platform  orientation  and  the  general 
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definition  is  given  by  [3] 


k  d 


T  = 
np 


np 


(3.78) 


for  the  stationary  model.  There  are  now  M  target  unit  vectors  and  M  platform 
orientations,  thereby  creating  M  differential  time  delays 


,  _  k-(0-,0").d- 

—ndx  cos  0"  ©  sin  0’^'"  —  pdz  sin 
c 


(3.79) 

(3.80) 


Analogously  to  the  stationary  model  development,  spatial  frequencies  are  defined  from 
the  differential  time  delay  as 

~  rv  dx  cos  0’'''  0  sin  0'”'' 


(3.81) 


and 


~  rv  dz  sin  0^ 
=  - — 


(3.82) 


Physically,  Equation  (3.81)  describes  the  a:'’''-axis  phase  behavior  for  each  platform 
orientation  while  Equation  (3.82)  describes  the  .2'^''-axis  phase  behavior  for  each  plat¬ 
form  orientation.  The  stationary  model  azimuth  and  elevation  steering  vectors  are 
given  by  the  A”  x  1  and  P  x  1  vectors  [3] 


and 


(3.83) 


e{^z)  = 


(3.84) 
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Similar  steering  vectors  can  accurately  model  phase  across  elements  in  the  maneuver 
model,  but  only  for  a  single  platform  orientation: 


(3.85) 


and 

1^)  =  ei27n?7U(0:P-i)^_  ^3_gg) 

The  formatting  difficulty  arises  because  there  are  now  M  unique  azimuth  steering 
vectors  like  a  and  M  unique  elevation  steering  vectors  like  e  needed  to  model  phase 
across  elements.  The  e  ®  b  ®  a  steering  vector  paradigm  is  therefore  abandoned. 
Instead,  the  M  spatial  frequencies  describing  the  and  ;2''''-axis  phase  behavior  are 
recast  into  MNP  x  1  spatial-temporal  steering  vectors.  These  MNP  x  1  spatial- 
temporal  steering  vectors  contain  the  phase  contribution  due  to  physical  element 
separation  for  each  element  and  pulse.  The  per  element  per  pulse  phase  contribution 
due  to  horizontal  element  is  created  by, 

(3.87) 

where  N  —  \  N  —  1]^  and  Ip  is  a  P  x  1  vector  of  ones.  And  the  per  element  per 

pulse  phase  contribution  due  to  vertical  element  separation  is  created  by 


g  (.^rv  )  gj27rP®i?pi8ilN 


(3.88) 


where  P  =  [0  :  P  —  1]^  and  In  is  a  x  1  vector  of  ones.  The  total  spatial-temporal 
phase  due  to  element  separation  is  simply  e  ©  a.  The  final  component  needed  to 
capture  the  total  spatial-temporal  phase  of  the  return  signal  is  the  Doppler  phase 
contribution. 


3.J^.2  Temporal- Spatial  Steering  Vector  with  Motion.  The  temporal  phase 
delay  depends  on  the  pulse  to  pulse  Doppler  frequency  of  the  returns,  just  as  it  did 
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in  the  stationary  model.  The  normalized  Doppler  frequency  uj  is  defined  [9], 


=  if .  (3.89) 

Ajrjr 

where  Vr  is  the  relative  velocity  between  the  aircraft  and  target,  A  is  the  radar  transmit 
wavelength,  and  fr  is  the  PRF.  The  relative  velocity  between  two  objects  is  the  time 
rate  of  change  of  the  distance  between  them. 


Vr  = 


dR 

dt 


(3.90) 


This  time  rate  of  change  is  the  projection  of  the  aircraft  velocity  vector  onto  the 
unit  vector  to  the  target,  k,  plus  the  projection  of  the  target  velocity  vector  Vj  onto 
k.  There  are  now  M  unit  vectors  to  the  target  and  M  aircraft  velocity  vectors  due 
to  platform  maneuver,  resulting  in  M  different  relative  velocities 


Vr- 


k^v  .  v7  +  k"  •  Vf 


(3.91) 


and  M  different  Doppler  frequencies 


=  2- 


A/p 


(3.92) 


In  the  stationary  model,  the  relative  velocity  between  a  location  (be  it  a  target  or 
clutter)  remained  constant  throughout  the  CPI.  The  Doppler  frequency,  u,  was  con¬ 
stant  and  the  resulting  pulse  to  pulse  phase  due  to  Doppler  was  a  linear  function  of 
a  single  Doppler  frequency: 


b(c9)  =  ■ 


(3.93) 


In  the  maneuver  model,  the  M  different  relative  velocities  lead  to  M  different  Doppler 
frequencies  and  the  pulse  to  pulse  phase  change  due  to  Doppler  is  no  longer  a  linear 
function  of  a  single  normalized  doppler  value.  The  M  Doppler  frequency  values  given 
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by  Equation  (3.92)  correspond  to  individual  Doppler  frequencies, 


d;  =  [^’■'^(0)  ^■'''(1)  . . .  d"(M  -  1)]'^. 


(3.94) 


The  corresponding  pulse  to  pulse  phase  progression  due  to  these  individual  Doppler 
frequencies  is  sought.  Capturing  that  phase  progression  requires  the  M  individual 
Doppler  values  recast  in  the  form, 

M-l 

W  =  [d"(0)  +  d;"(l)]  ...  (3.95) 

m=0 


Notice  Equation  (3.94)  collapses  to  the  stationary  model  linear  form  [cu  2u  ... 
if  a;(0)  =  C(;(l)  =  A  convenient  matrix  algebra  structure  is  created  to 

transform  the  M  Doppler  values  in  Equation  (3.94)  into  the  pulse  to  pulse  Doppler 
progression  of  Equation  (3.95).  An  upper  triangular  matrix  of  u)  values  is  first  created. 


cu"(0)  cu"(l)  d"(2) 


d"(M  -  1) 


0  d"(0)  d"(l) 

0  0  ^’^'^(0) 

0  0 

0  0 


d"(M  -  2) 


cu'''^(0) 


where  the  subscript  represents  the  pulse  number  within  the  CPI.  Multiplying  O  by  a 
ones  column  vector  of  length  M  yields  the  progression  given  in  Equation  (3.95), 


=  DIm  (3.96) 

M-l 

=  [d"(0)  [d"(0)+(h'^^(l)]  ...  5]]  d'^^(m)]^.  (3.97) 

m=0 

And  the  corresponding  pulse  to  pulse  Doppler  phase  contribution  is  the  M  x  1  vector 
gi27rw  Although  Equation  (3.97)  produces  a  M  x  1  column  vector  dimensionally  the 
same  as  the  stationary  model,  the  M  x  1  vector  must  be  recast  into  a  MNP  vector 
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to  match  and  0(1?").  The  per  element  per  pulse  phase  contribution  due  to 

Doppler  element  is  created  by, 


b(w)  =  ('3_gg) 

And  the  complete  steering  vector  with  motion  to  location  (d,0)  is  therefore  given  by 

v  =  e(^r)©b(w)©ara,  (3.99) 

where  e  represents  the  spatial-temporal  phase  contribution  due  to  vertical  element 
separation,  b  represents  the  temporal-spatial  phase  contribution  due  to  Doppler,  and 
a  represents  the  spatial-temporal  phase  contribution  due  to  horizontal  element  sepa¬ 
ration.  Both  the  maneuver  steering  vector  of  Equation  (3.99)  and  stationary  steering 
vector  of  Equation  (2.26)  are  MNP  x  1. 

The  introduction  of  platform  maneuver  resulted  in  M  values  for  previously  de¬ 
fined  single  value  characteristics.  New  MNPxl  spatial-temporal  (e,  a)  and  temporal- 
spatial  (b)  steering  vectors  replaced  the  original  spatial  (e,  a)  and  temporal  (b) 
steering  vectors  of  the  stationary  model.  The  e  ©  b  ©  a  steering  vector  paradigm  was 
subsequently  replace  by  e  0  b  ©  a.  The  final  piece  needed  for  a  space-time  snapshot 
remains  the  return  amplitude  with  motion. 

3.4-3  Return  Amplitude  with  Motion.  The  complete  space-time  snapshot  to 
a  location  is  the  return  amplitude  times  the  steering  vector, 

X  =  ttv.  (3.100) 

Previously  a  was  assumed  constant  throughout  the  CPI.  The  assumption  was  viable 
because  the  environment  was  identically  illuminated  throughout  the  CPI.  Platform 
maneuver  invalidates  that  assumption  as  locations  are  now  illuminated  differently  as 
the  platform  rotates.  The  amplitude  must  be  recalculated  each  pulse,  resulting  in  a 
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M  X  1  vector  of  amplitudes.  The  actual  calculations  of  these  amplitudes  are  discussed 
in  the  target  model  (Section  3.5)  and  clutter  model  (Section  3.6).  The  formatting  of 
these  M  X  1  vector  of  amplitudes  to  match  the  maneuver  steering  vector  framework 
is  discussed  here.  With  the  M  x  1  vector  of  amplitudes  given  by  o:,  the  MNP  x  1 
amplitude  vector  is 

d  =  Ip  (g)  a  O  In-  (3.101) 

And  the  complete  space-time  snapshot  with  platform  maneuver  is  given  by 

X  =  d  ©  e(i?7)  ©  b(w)  O  a(i?^'').  (3.102) 

The  target  and  clutter  models  are  next  developed  incorporating  platform  maneuver. 
Note  the  four  M  x  1  components  required  to  build  the  target  and  clutter  models:  re¬ 
turn  amplitude  a,  spatial  frequencies  and  ^  and  normalized  Doppler  frequency 
u}.  Once  these  M  x  1  quantities  are  characterized,  the  formatting  methodology  built 
in  this  section  produces  the  space-time  snapshot.  Equation  (3.101)  formats  the  M 
return  amplitudes  o:  into  the  MNP  x  1  d.  Equation  (3.87)  formats  the  M  hori¬ 
zontal  spatial  frequencies  into  the  MNP  x  1  spatial-temporal  steering  vector  a(i9"). 
Equation  (3.88)  formats  the  M  vertical  spatial  frequencies  into  the  MNP  x  1  spatial- 
temporal  steering  vector  e('i?^'').  Equations  (3.94-3.97)  format  the  M  normalized 
Doppler  frequencies  into  the  MNP  x  1  temporal- spatial  steering  vector  b(w).  And 
finally  Equation  (3.102)  combines  the  three  MNP  x  1  steering  vectors  capturing  the 
return  phase  and  one  MNP  x  1  steering  vector  capturing  the  return  magnitudes  into 
the  complete  space-time  snapshot. 

3. 5  Target  Model 

The  four  M  x  1  components  required  to  build  the  target  model  are:  return 
amplitude  ott,  spatial  frequencies  i?"  and  and  normalized  Doppler  frequency  ujf 
The  target  location  is  defined  as  {Ot,  (t>t)  as  measured  from  the  stationary  platform 
orientation.  Platform  maneuver  yields  M  azimuth  and  elevation  angles  on  transmit, 
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and  M  on  receiver,  {6’^ ,4)^).  The  stationary  model  Signal-to-Noise- Ratio 
(SNR)  is  defined  as  [3] 

,  _  PtG{0t,4t)g{et,4t)X^at 

(AwYNBLsRtiOy  ’  ^  ^ 

where  Pt  is  transmit  power,  G  is  transmit  array  gain,  g  is  element  pattern  power  on 
receive,  A  is  transmit  wavelength,  at  is  target  radar  cross  section  (RCS),  NBLg  are 
noise  and  system  losses,  and  Rt  is  range  to  target.  In  the  stationary  model’s  imple¬ 
mentation  the  SNR  is  set  to  unity.  In  both  models  the  target  snapshot  amplitudes 
obey  [3] 

£{atat}  =  (3.104) 


where  is  the  noise  power.  A  constant  SNR  in  the  stationary  model  is  valid  be¬ 
cause  the  target  is  illuminated  identically  throughout  the  CPI.  However,  a  constant 
maneuver  model  SNR  fails  to  capture  the  true  target  return  magnitude.  A  tar¬ 
get  may  reside  in  the  mainbeam,  sidelobes,  and  nulls  within  a  single  CPI  and  re¬ 
turns  from  pulse  to  pulse  may  have  varying  SNR.  Returning  to  the  stationary  model 
SNR  definition  of  Equation  (3.103),  the  only  platform  maneuver  induced  changes 
are  the  transmit  array  gain  G{9t,(t>t)  and  receive  element  pattern  g{0t,(f>t),  introduc¬ 
ing  M  values  of  each  quantity.  Compared  to  the  stationary  model,  the  maneuver 
model  SNR  will  be  greater  if  G{6^^ ,(j)^^)g{6'^^ ,4)^)  >  G{9t,4>t)g{dt,4t),  and  less  if 
G{0^^,  4‘^^)g{0^'' ,  4R')  <  G{9t,  4>t)g{9t,  4t)-  The  maneuver  model  SNR  is  therefore  de¬ 
fined  as  a  scaled  function  of  the  stationary  model  G{9t,  4>t)g{(^t,4t), 


•  G(0^-,</)^-)0g(0-,O 

G{9t,  (j)t)g{(^t^  4t) 


(3.105) 


Anytime  G{0^^,4>^^)g{0^^,4>^^)  =  G{9t,4>t)g{dt,4t),  the  target  is  illuminated  exactly 
the  same  in  both  models  and  the  maneuver  SNR  will  be  equal  to  the  stationary 
SNR.  Anytime  G{0^^,  4>*'^)g{0^^ ,  4>'^^)  >  G{9t,4>t)g{9t,4’t),  the  target  is  illuminated  at 
a  higher  power  in  the  maneuver  model  and  the  SNR  will  be  greater  than  the  stationary 
SNR.  And  anytime  G{0^^,  0*’^)p(0",  4>^^)  <  G{9t,  4>t)g{0t:  4t)i  the  target  is  illuminated 
at  a  lower  power  in  the  maneuver  model  and  the  SNR  will  be  less  than  the  stationary 
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SNR.  The  M  x  1  SNR  values  of  Equation  (3.105)  are  reshaped  into  the  MNP  x  1 
values  by 

=  (3.106) 


And  by  the  dehnition  in  Equation  (3.104),  the  MNP  x  1  vector  of  target  amplitudes 
is 


oct  =  a 


(3.107) 


With  the  target  amplitude  characterized  in  Equation  (3.107),  the  next  target  model 
component  requiring  characterization  is  normalized  Doppler.  Target  velocity  may 
include  x,  y,  and  z  components,  denoted  Vt^,  and  Vt^. 


Vt  =  Vt,x  +  Vt^y  +  VtJ 


(3.108) 


The  normalized  Doppler  frequency  for  a  target  at  was  present  in  Equation 

(2.45)  as 


Ut  —  2 


(va  +  v^)  cos  9t  sin  (pt  +  vf  cos  6t  cos  pt  +  sin  Ot 


(3.109) 


No  consideration  of  the  target  direction  of  travel  was  necessary  in  the  stationary 
model.  Because  the  relative  velocity  between  the  aircraft  and  target  was  constant, 
the  Doppler  frequency  was  constant  and  u  could  simply  be  assigned  a  value  that  would 
not  change  during  the  CPI.  When  platform  maneuver  is  induced,  a  target  will  possess  a 
different  Doppler  frequency  at  each  platform  orientation.  Target  velocity  is  restricted 
to  VtyY  in  this  model’s  implementation,  corresponding  to  the  target  travelling  toward 
the  aircraft  along  the  y-axis  at  speed  Vty  and  given  by 


Vt  =  Vtyj. 


(3.110) 


The  normalized  Doppler  expression  with  restricted  target  velocity  becomes 


LVf  —  2 


Va  cos  9t  sin  pt  +  v\  cos  9t  cos  pt 


(3.111) 
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In  the  production  of  SINK  plots,  a  target  normalized  Doppler  is  specified.  The  ma¬ 
neuver  target  model  will  therefore  calculate  the  target  velocity  required  to  produce 
the  original  normalized  Doppler  value  of  interest  by 


OJtXfr  -  Va  COS  Of  siu 


COS  6t  cos  (/)f 


(3.112) 


And  the  M  normalized  Doppler  values  for  the  target  at  then  become 


—  2 


Va  COS  0’'''  ©  sin  (j/^  +  vl  cos  0  COS0’' 
^fr 


(3.113) 


The  target  temporal  steering  vector  with  motion  is  constructed  with  the  upper  trian¬ 
gular  matrix 

^  Ot{2)  ...  OtiM-1) 

0  cUi(O)  (hi(l)  ...  LDt{M-2) 

0  0  u!tio)  ...  ; 

0  0  ...  •••  ; 

0  0  .  h;i(0) 

multiplied  by  the  M  x  1  vector  of  ones 


Wj  —  Sljljvr, 


(3.114) 


and  recast  into  the  MNP  x  1  snapshot  sequence  by 


(3.115) 


The  final  two  pieces  of  the  target  space-time  snapshot  are  the  spatial-temporal 
steering  vectors  describing  the  phase  contribution  due  to  element  separation.  The 
M  X  1  spatial  frequencies  are  given  by 

~  rv  dx  COS  0’'''  ©  sin 


(3.116) 
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and 


~  rv  dz  sin  O' 
'd  =  — - 

c  ’ 


(3.117) 


where  0"  and  are  the  M  radar  angles  to  the  target  at  (d,0).  The  corresponding 
steering  vectors  are  given  by 


rv\  _  j27rlp(8)'i?^'^(g)N 


(3.118) 


and 


rv\  _  j27rP(g)i?pi8ilisr 


=  e 


(3.119) 


And  the  complete  target  snapshot  can  now  be  assembled  as 


Xt  =  at(D  et(^7)  ©  bt(wt)  0  at(i?") 


(3.120) 


3.6  Clutter  Model 


Accurate  characterization  of  the  clutter  environment  requires  a  return  amplitude 
otik,  spatial  frequencies  6^  and  6^,  and  normalized  Doppler  frequency  uhfc  for  each 
ik  patch.  The  clutter  location  is  defined  as  (9ik,  (pik)  as  measured  from  the  stationary 
platform  orientation.  Platform  maneuver  yields  M  azimuth  and  elevation  angles  on 
transmit,  (0*^,0^^),  and  M  on  receiver, 


The  patch  is  now  illuminated  differently  each  pulse,  requiring  per  pulse  calcula¬ 
tion  of  Clutter-to-Noise- Ration  (CNR)  ^ik  and  amplitude  ctj^.  The  stationary  model 
CNR  is  defined  as  [3] 


PtG{et,(Pt)9{0t,<Pt)X^aik 

{ATTfNBLMikiOy 


(3.121) 


where  Pt  is  transmit  power,  G  is  transmit  array  gain,  g  is  element  pattern  power 
on  receive,  A  is  transmit  wavelength,  cTj^  is  clutter  patch  RCS,  NBLg  are  noise  and 
system  losses,  and  Rik  is  range  to  the  patch.  The  patch  RCS  definition  of  Equation 
(2.63)  remains  unchanged  because  the  patch  area,  grazing  angle,  and  reflectivity  are 
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unchanged.  The  M  CNR  values  are  given  by 


^ik  ^ik 


{AnfNBLsRikiey 


and  are  recast  into  MNP  x  1  column  vector  by 


lifc  —  Ip  ®  kik  ®  In- 


(3.122) 


(3.123) 


The  complex  amplitude  6;,^  obeys 


^{^ikm^jln} 


=  a 


^ikm^jln^i—j^k—l  i 


(3.124) 


meaning  the  returns  from  different  clutter  patches  are  uncorrelated.  Next  consider 
the  M  Doppler  values  a  single  patch  and  the  construction  of  the  M NP  x  1  temporal- 
spatial  steering  vector  that  captures  the  phase  contribution  due  to  Doppler. 


A  clutter  patch  was  previously  located  at  one  d,  one  (f>,  one  u,  and  one  return 
amplitude  ac-  The  normalized  Doppler  frequency  for  a  patch  at  {0,(1))  was  presented 
in  Equation  (2.55)  as 


UJik  2 


cos  0i]„  sin 

Wr 


(3.125) 


The  M  normalized  Doppler  values  for  a  target  at  (6*,0)  are  given  by 


—  2 


Va  cos  0^^  ©  sin 

Wr 


(3.126) 
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The  clutter  patch  temporal  steering  vector  with  motion  is  constructed  with  the  upper 
triangular  matrix 


^ik  — 


^ifc(O)  ^ifc(2) 

0  ^i/c(0)  (Ujfc(l) 

0  0  iVik{0) 

0  0 

0  0 


LUik{M  -  1) 
^ik{M  —  2) 


multiplied  by  the  M  x  1  vector  of  ones 


M- 


And  recast  into  the  MNP  x  1  form  by 


(3.127) 


(3.128) 


The  hnal  pieces  needed  for  the  clutter  space-time  snapshot  are  the  spatial- 
temporal  steering  vectors,  given  by 

^3_^29) 


and 


(3.130) 


The  antenna  array  radiates  in  all  directions,  illuminating  every  i  x  k  clutter 
patch.  Therefore,  the  complete  clutter  space-time  snapshot  is  the  sum  oiixk  clutter 
patch  snapshot 

Nr-l  Wc-1 

'^^ikQ  e(i97)  ©  b(wifc)  0  ai-dYk)  (3.131) 

i=0  k=0 
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And  the  clutter  covariance  matrix  with  platform  maneuver  is  given  by, 


Rc  =  ^{XcXc) 


(3.132) 


Nr-l  Nc-1  / 


(3.133) 


3.7  Maneuver  Model  Summary 

This  Chapter  developed  a  maneuver  model  based  on  the  popular  stationary 
model  of  [3, 10].  The  maneuver  model  developed  can  simulate  any  degree  of  pitch, 
roll,  yaw,  or  any  combination  of  the  three.  This  research  focuses  on  development  of  the 
model  itself,  and  therefore  the  model  investigates  individual  rotations  versus  specific 
innovative  and  complex  maneuver  scenarios.  The  overriding  theme  of  the  maneuver 
model  is  that  singular  quantities  in  the  stationary  model  of  [3,10]  possess  M  unique 
quantities  in  the  maneuver  model.  The  platform  defines  the  reference  coordinate  sys¬ 
tem  at  all  times.  The  {6,  <p)  radar  coordinates  of  each  clutter  patch  and  the  target  are 
first  transformed  into  their  cartesian  counterparts  {x,y,z).  Rotation  matrices  are  used 
to  recalculate  the  coordinates  to  each  {x,y,z)  point  based  on  the  amount  of  rotation, 
measured  with  the  platform  orientation  as  the  reference.  The  resulting  {x',y',z')  coor¬ 
dinates  are  transformed  back  to  radar  coordinates  (9',  <p').  Each  location  has  M  {O' , 
(f)')  pairs  on  transmit  and  receive,  translating  to  M  Doppler  frequencies,  spatial  fre¬ 
quencies,  and  amplitudes.  Four  MNP  x  1  vectors  are  created,  one  for  amplitude,  one 
each  for  phase  contributions  due  horizontal  spacing,  vertical  spacing,  and  Doppler. 
The  Kronecker  product  paradigm  is  replaced  with  the  Hadamard  product  in  creating 
the  space-time  snapshots.  The  maneuver  models  allows  characterization  of  the  target 
and  interference  (thermal  noise  and  clutter  for  this  research)  environment.  The  next 
logical  step  is  investigating  maneuver  impacts  on  the  interference  environment,  ME 
output  SINK,  and  sample  support  homogeneity.  These  topics  are  addressed  next  in 
Chapter  IV. 
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IV.  Maneuver  Effects  on  Clutter  Spectrum,  SINK,  and 

Homogeneity 


Using  the  model  developed  in  Chapter  III,  platform  maneuver  effects  are  simu¬ 
lated  for  the  side-looking  (SL)  airborne  radar.  The  objective  is  characterization 
of  maneuver  effects  on  clutter  spectrum,  Matched  Filter  (MF)  Signal-to-Interference- 
plus-Noise- Ratio  (SINK),  and  sample  support  homogeneity.  The  number  of  azimuth 
channels,  elevation  channels,  and  target  ranges  are  varied  to  gain  insight  into  maneu¬ 
ver  effects.  A  rotation  rate  of  800°/sec  and  pulse  length  of  32  pulses  is  used  unless 
otherwise  noted.  The  rotation  rate  corresponds  to  12.86°  total  rotation  during  the 
CPI,  similar  to  the  forward-looking  (FL)  work  in  [4].  The  number  of  azimuth  and 
elevation  channels  is  varied  on  a  case  by  case  basis  to  illustrate  maneuver  effects. 
The  focus  is  airborne  moving  target  indication,  thus  the  array  is  steered  to  boresight 
(^  =  0,  0  =  0).  The  array  is  steered  at  boresight  throughout  the  coherent  processing 
interval,  but  its  orientation  relative  to  the  clutter  and  target  changes  as  the  platform 
rotates.  Two  key  phenomenon  identified  are  strongly  correlated  with  Matched  Fil¬ 
ter  (MF)  Signal-to-Interference-plus-Noise  (SINK)  and  sample  support  homogeneity. 
The  first  is  the  maneuver  induced  differences  in  clutter  illumination,  and  the  second 
is  the  maneuver  induced  differences  in  target  illumination.  As  each  maneuver  type 
is  examined,  visual  plots  depicting  these  (target  and  clutter  illumination)  differences 
are  used  to  illustrate  maneuver  effect  and  gain  insight.  Azimuthal  antenna  pattern 
cuts  illustrate  the  clutter  illumination  magnitudes  at  specific  ranges.  Signal  Match 
(SM)  power  spectral  density  (PSD)  and  minimum  variance  estimate  (MVE)  plots  il¬ 
lustrate  the  changes  in  the  clutter  environment.  Signal-to-Noise- Ratio  (SNR)  plots 
depict  changes  in  target  illumination.  These  characterizations  of  the  changing  inter¬ 
ference  environment  and  target  illumination  lead  directly  into  the  MF  SINR  and  data 
homogeneity  results.  The  results  are  organized  as  follows: 

Section  4.1  gives  a  brief  review  of  the  clutter  spectral  estimations  and  perfor¬ 
mance  metrics,  and  their  specific  implementation  within  the  stationary  and  maneuver 
models.  Section  4. 2-4.4  describes  the  pitch,  roll,  and  yaw  effects  for  linear  and  planar 
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arrays  at  two  target  ranges.  And  Section  4.5  describes  the  maneuver  induced  hetero¬ 
geneities  in  the  sample  support  vectors.  The  simulation  parameters  are  described  in 
Table  4.1. 


4-1  Spectral  Estimation  and  Performance  Metrics 

The  SM  spectral  estimation  method  was  presented  in  Chapter  II  and  is  applied 
here.  The  stationary  PSD  is  given  by 


P  = 


v^v 


(4.1) 


where 


V  =  0(^2)  ®  b((A;)  ®  a((A;), 

(4.2) 

R-  =  R-n  +  Rc) 

(4.3) 

and 

Re  =  e{xx^} 

Nr-l  Nc-l 

=  (7^  ®  b(n;)b(n;)^  ®  . 

i=0  k=0 

(4.4) 

The  maneuver  PSD  is  given  by 

~  v^Rv 

if  ’ 

V 

(4.5) 

where, 

R  =  +  Rc, 

(4.6) 
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Table  4.1:  Simulations  parameters. 


Variable 

Value 

M 

32 

N 

varied 

P 

varied 

fo 

1240  MHz 

fr 

1984  Hz 

Tp 

0.8  ps 

Pt 

200  kW 

B 

800  kHz 

Fn  (Noise  Figure) 

3  dB 

Nc 

361 

ha  (aircraft  altitude) 

3073  m 

Va  (aircraft  velocity) 

dxfr 

2 

R 

11,66  km 

7 

-3  dB 

Array  Transmit  Gain 

22  dB 

Element  Pattern 

Gosine 

Element  Gain 

4  dB 

Element  Backlobe  Level 

-30  dB 

dx 

0.10922  m 

dz 

0.1407  m 

Transmit  Taper 

Uniform  (None) 

System  Losses  Lg 

3  dB 

Target  (f) 

0° 

Target  9 

0° 

Be 

10 

fr 

Ge 

4  dB 

p  (rotation  rate) 

800° /sec 

(stationary) 

1 
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and 


Rc 


^{XcXc  } 

Nr-1  Nc-1 


a 


i=0  k=0 


Y  (y^ik®  ^i'^7)  ®  b(wifc)  © 


'hk  ©  e(i97)  0  b(wifc)  0  a(i?") 


H 


(47) 


(4^8) 


The  stationary  Minimum  Variance  Estimate,  discussed  in  Chapter  II,  for  the  station¬ 
ary  scenario  is 

P^ve  =  (v'^R-V)-^  (4.9) 

And  for  the  maneuver  model  is  given  by 


^mve  =  (V^R  'v) 


(4.10) 


The  stationary  MF  output  SINK  is  given  by 


|2 


SINR  = 


a^t  I  w"vf 

w^Rw  ’ 


(4.11) 


where  w  is  the  non-motion  compensated  steering  vector  defined  as 


w  =  R“^v, 


(4.12) 


and  is  set  to  unity.  And  the  maneuver  MF  output  SINK  is  given  by 


SINR”^  = 


I  w^(^t  0  vt) 

w^Rw 


(4.13) 


where 


w  =  R"V 


(4.14) 
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-2.74 


Figure  4.1:  A  12.86°  pitch  rotation  results  in  only  —0.63°  change  in  azimuth  angle 
and  +0.07°  change  in  elevation  angle  to  the  location  0  =  0°  at  a  range  of  66  km. 

4 '2  Pitch  Effects 

Platform  pitch  corresponds  to  rotation  about  the  y-axis.  The  antenna  is  steered 
to  boresight,  where  the  target  is  located  throughout  the  CPI.  In  the  pitch  scenario, 
the  target  is  always  located  at  boresight  despite  the  platform  maneuver.  The  location 
of  the  y-axis,  defining  boresight  and  target  location,  never  changes  as  pitch  refers  to 
rotation  about  that  axis.  Therefore,  the  target  resides  in  the  main  beam  at  boresight 
at  all  times,  and  is  constantly  illuminated  =  1m). 

4-2.1  Linear  Array  with  Target  at  11  km  and  66  km.  Pitch  changes  both  the 
azimuth  and  elevation  angles  to  the  location  0  =  0  at  66  km,  as  depicted  in  Figure 
4.1.  Notice  a  12.86°  pitch  rotation  corresponds  to  changes  in  azimuth  and  elevation  of 
just  —0.63°  and  +0.07°  degrees  respectively.  Significant  SINK  and  detection  impact 
are  therefore  expected  only  if  the  array  pattern’s  azimuth  and  elevation  definition 
significantly  changes  over  that  short  span.  In  the  case  of  the  linear  array,  the  pattern 
lacks  elevation  definition,  and  a  3  dB  azimuthal  beamwidth  on  the  order  of  0.63° 
is  not  achieved  unless  N  is  at  least  115  elements  long.  Figure  4.2  illustrates  an 
antenna  pattern  azimuth  cut  at  66  km  for  a  120  x  1  array  for  the  stationary  array 
orientation  and  for  12.86°  pitch  orientation,  corresponding  to  the  antenna  pattern 
on  the  final  pulse  of  the  CPI.  Shorter  linear  arrays  lack  the  azimuthal  definition  to 
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Figure  4.2:  The  small  change  in  azimuth  and  elevation  angles  at  66  km  due  to  pitch 
requires  greater  definition  in  the  antenna  pattern  to  illuminate  the  ground  with  any 
significant  difference.  A  120  x  1  array  pattern  is  projected  across  azimuth  at  66  km, 
demonstrating  the  need  for  many  elements  to  achieve  any  discernable  difference. 

create  significant  changes  in  the  clutter  illumination  at  66  km.  Consequently,  the 
clutter  environment  remains  relatively  unchanged  until  the  array  length  approaches 
115  X  1.  Figure  4.3  contains  four  plots  depicting  SM  PSDs  and  MVEs  for  both  the 
stationary  and  pitch  scenarios  at  66  km.  Note,  the  MVE  plots  provided  are  for  a 
64  X  1  linear  array.  The  azimuthal  and  Doppler  resolution  needed  to  create  complete 
PSDs  for  a  120  x  1  array  is  processing  prohibitive  and  differences  are  nearly  impossible 
to  see.  The  64  x  1  result  provides  ample  physical  insight  into  the  scenario  and  can  be 
extended  to  the  120  x  1  array.  Overall,  pitch  did  not  have  an  overwhelming  effect  on 
the  clutter  spectral  content.  It  appears  the  clutter  ridge,  the  “s” -shaped  ridge  from 
negative  azimuth  and  Doppler  to  positive  azimuth  and  Doppler,  shifted  slightly  in 
the  negative  Doppler  direction.  An  azimuthal  SM  cut  at  0  =  0  shown  in  Eigure  4.4 
confirms  the  slight  negative  Doppler  shift  in  the  clutter  spectrum.  The  SM  PSD  cut 
in  Eigure  4.4  directly  translates  to  the  output  SINK  plot  in  Eigure  4.5.  The  increased 
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Figure  4.3:  Spectral  estimation,  64  x  1  array  with  pitch  at  66  km.  MVE  for 

the  stationary  (top  left)  and  pitch  (top  right)  64  x  1  arrays  shows  slight  spread  in 
the  clutter  spectrum  induced  by  the  pitch  rotation.  A  slight  shift  towards  negative 
Doppler  is  noted  in  the  SM  PSD  pitch  scenario  (bottom  right),  compared  to  the  SM 
PSD  stationary  (bottom  left).  This  shift  is  later  confirmed  with  a  SM  PSD  cut  at  0 
azimuth. 

negative  clutter  power  around  zero  Doppler  translates  to  the  SINK  degradation  when 
pitch  is  introduced,  labeled  with  a  dashed  line  and  diamond  (o)  in  Figure  4.5.  Notice 
the  improvement  between  MF  and  NA  performance  increases  significantly.  The  MF 
improvement  over  NA  at  0.1  normalized  Doppler  was  ~  20  dB  in  the  stationary 
scenario.  That  improvement  over  NA  increased  to  ~  30  dB  with  pitch  maneuver, 
demonstrating  that  NA  processing  is  much  more  susceptible  to  pitch  than  the  MF. 

If  the  target  is  located  at  11  km  versus  66  km,  pitch  effects  are  more  pronounced. 
Figure  4.6  shows  the  azimuth  and  elevation  angles  to  0  =  0  changes  more  rapidly  as 
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Figure  4.4:  SM  PSD  cut  at  zero  azimuth,  120  x  1  with  pitch.  A  slight  change  in  the 
SM  PSD  is  noted  for  pitch  at  66  km.  The  shift  toward  negative  Doppler  translates 
to  the  null  shift  in  the  output  SINK  plot. 

range  decreases.  A  12.86°  pitch  rotation  for  0  =  0  at  11  km  corresponds  to  changes 
in  azimuth  and  elevation  of  —3.35°  and  +0.04°  respectively,  a  considerable  difference 
from  the  —0.63°/  +  0.07°  at  66  km.  The  number  of  azimuth  elements  required  to 
impact  the  clutter  environment  and  SINK  is  reduced  from  115  at  66  km  to  25  at 
11  km.  In  order  to  compare  results  between  11  and  66  km,  the  120  x  1  array  is  used 
again.  The  SM  PSD  cut  for  0  =  0  at  11  km  is  shown  in  Figure  4.8.  The  resulting 
spectral  estimates  in  Figure  4.7  show  the  dramatic  change  in  the  clutter  picture. 
Compared  with  the  cut  for  0  =  0  at  66  km  (Figure  4.4),  the  negative  Doppler  shift  is 
much  more  pronounced  and  all  Doppler  lobing  structure  is  lost.  The  mainbeam  width 
for  the  120  x  1  array  is  roughly  0.06°  (see  Figure  4.2).  An  azimuth  sweep  of  3.35°, 
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Figure  4.5:  Output  SINK  for  120  x  1  array  with  pitch  at  66  km.  The  shift  in  SM 
PSD  towards  negative  Doppler  is  seen  in  the  output  SINK.  Overall,  the  array  pattern 
azimuth  definition  is  just  approaching  the  definition  needed  to  cause  degradation. 
Improvement  over  NA  processing  with  pitch  increased  considerably. 

corresponding  to  the  change  in  azimuth  at  11  km,  sweeps  through  nearly  three  lobes 
of  the  antenna  pattern.  The  clutter  illumination  for  a  location  at  11  km  experiences 
much  greater  fluctuation  than  one  at  66  km,  and  the  resulting  Doppler  spread  seen  in 
Figure  4.8  results.  The  changing  clntter  environment  depicted  in  the  PSD  translates 
to  the  output  SINK  in  Figure  4.9.  First,  the  widened  negative  clutter  null  in  the  PSD 
is  present  in  the  NA  processing  SINK.  Next,  it  is  also  observed  in  the  widening  of  the 
main  clutter  null  in  the  MF.  The  clntter  nnll  actually  shifted  to  negative  Doppler, 
away  from  zero  Hz.  Similar  to  the  resnlt  at  66km,  the  MF  improvement  gap  over  NA 
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Figure  4.6:  A  12.86°  pitch  rotation  results  in  a  —3.35°  change  in  azimuth  angle  and 
+0.04°  change  in  elevation  angle  to  the  location  0  =  0°  at  a  range  of  11  km.  The 
change  in  azimuth  at  11  km  is  nearly  3  full  degrees  more  than  that  at  66  km. 

processing  widened  when  pitch  was  introduced.  The  increased  clutter  power  due  to 
smaller  range  resulted  in  MF  improvement  of  ~  45  dB  over  NA  processing. 

Next,  pitch  effects  are  considered  for  an  increasing  number  of  azimuth  channels. 
As  the  number  of  azimuth  channels  are  increased,  the  transmitted  antenna  pattern 
gains  azimuthal  definition.  If  the  azimuth  resolution  is  high  enough,  the  resulting 
clutter  spectrum  broadens  and  shifts.  If  the  resolution  is  not  high  enough,  the  clutter 
spectrum  remains  close  to  the  stationary  spectrum.  Figure  4.10  shows  SM  cuts  at 
zero  azimuth  for  N  =  100,  150,  and  200.  Note,  the  CPI  was  shortened  to  20  pulses 
for  these  plots  to  computationally  allow  increased  N  values.  Severe  spreading  is  seen 
for  the  11km  scenario,  while  the  66  km  case  is  relatively  stable  due  to  the  small 
changes  in  azimuth  and  elevation  angles  pitch  causes  at  that  range.  If  the  azimuth 
resolution  is  too  low  for  significant  pitch  effects,  increasing  azimuth  channels  improves 
MF  performance.  This  result  is  observed  in  the  narrowing  null  in  the  SINK  loss  plot 
of  Figure  4.11.  On  the  other  hand,  if  the  azimuth  resolution  is  high  enough  such  that 
pitch  effects  are  significant,  increasing  azimuth  channels  results  in  the  clutter  null 
spread  and  shift  seen  in  Figure  4.11. 

Planar  Array  with  Target  at  11  km  and  66  km.  Adding  elevation 
channel  changes  the  antenna  pattern  and  the  clutter  illumination.  The  entire  array 
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Figure  4.7:  Spectral  estimation,  64  x  1  array  with  pitch  at  11  km.  MVE  for 

the  stationary  (top  left)  and  pitch  (top  right)  64  x  1  arrays  shows  slight  spread  in 
the  clutter  spectrum  induced  by  the  pitch  rotation.  A  slight  shift  towards  negative 
Doppler  is  noted  in  the  SM  PSD  pitch  scenario  (bottom  right),  compared  to  the  SM 
PSD  stationary  (bottom  left).  This  shift  is  later  confirmed  with  a  SM  PSD  cut  at  0 
azimuth. 


gain  is  normalized  to  22  dB,  regardless  of  the  total  number  of  elements.  Therefore 
the  target  (at  boresight)  remains  under  the  same  illumination  as  the  linear  array. 
However,  the  clutter  is  illuminated  with  less  power  by  elevation  sidelobes  previously 
not  present  in  the  linear  array.  The  difference  in  illumination  is  more  pronounced 
in  the  11  km  scenario  because  of  the  elevation  angle  to  the  range.  Referring  back 
to  Figures  4.1  and  4.6,  the  elevation  angle  to  0  =  0  at  66km  is  —2.8°  and  to  11km 
is  —14°.  Figure  4.12  illustrates  this  concept  that  the  degree  of  impact  on  clutter 
illumination  depends  on  the  range.  The  addition  of  five  (5)  elevation  channels  results 
in  a  12  dB  reduction  in  clutter  illumination  at  11km,  but  only  a  0.3  dB  reduction 
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Figure  4.8:  SM  PSD  cut  for  120  x  1  array  at  0  =  0  for  11  km  with  pitch.  The 

increased  azimuth  resolution  of  the  antenna  pattern  and  the  change  in  azimuth  asso¬ 
ciated  with  pitch  at  11  km  create  the  spreading  and  shift  observed. 

at  66km.  The  illustration  at  66km  is  omitted  for  the  insignificant  change.  The 
reduction  in  clutter  illumination  at  1 1km  translates  to  a  reduction  of  the  clutter  power 
spectrum,  as  depicted  in  Figure  4.13.  Not  surprisingly,  the  clutter  power’s  significant 
reduction  in  the  planar  array  plus  the  added  DOF  translate  to  improved  SINK,  shown 
in  Figure  4.15.  While  increasing  azimuth  channels  increased  pitch  effects,  increasing 
elevation  channels  decreases  pitch  effects.  The  added  elevation  definition  reduces 
clutter  illumination  and  therefore  power,  while  the  target  is  illuminated  constantly  in 
the  normalized  mainbeam  at  22  dB. 

A  brief  recap  is  needed  of  pitch  effects  for  linear  and  planar  arrays,  including 
targets  at  close  (11km)  and  distant  (66km)  ranges.  Pitch  effects  are  insignificant  at 


Figure  4.9:  Output  SINK  for  120  x  1  array  with  pitch  at  11  km.  The  broadening 
and  larger  shift  in  SM  PSD  towards  negative  Doppler  is  seen  in  the  output  SINK. 
Improvement  over  NA  processing  again  increased,  just  as  it  did  at  66km. 

distant  ranges  for  the  linear  array,  as  long  as  the  array  length  and  resnlting  azimnthal 
definition  does  not  approach  the  corresponding  change  in  azimnth  resulting  from 
the  pitch  rotation.  For  the  simulation  parameters  used,  linear  array  pitch  effects  at 
66km  were  minimal  until  the  array  length  approached  115  elements.  The  greater 
the  target  distance,  the  greater  the  nnmber  of  elements  required  before  pitch  causes 
SINK  degradation.  Before  the  number  of  azimuth  channels  reaches  that  threshold 
number,  adding  azimuth  channels  improves  performance  as  if  pitch  were  not  present. 
Adding  elevation  channels  offers  marginal  improvement  for  distant  targets.  The  small 
elevation  angles  to  distant  targets,  and  the  very  small  change  in  elevation  angle  pitch 
induces,  requires  the  addition  of  many  elevation  channels  to  impact  the  results.  For 
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Normalized  Doppler 


Figure  4.10:  SM  PSD  cut  for  linear  arrays  at  0  =  0  for  11  km  and  66  km  with  pitch. 
The  number  of  pulses  in  the  CPI  was  reduced  to  20,  in  order  to  avoid  computation 
memory  limits.  Top  left:  SM  cut  for  11  km  without  platform  maneuver.  Top  right: 
SM  cut  for  11  km  with  pitch.  As  the  azimuth  channels  are  increased,  the  resolntion 
of  the  antenna  pattern  improves.  These  improvements  in  resolntion  lead  to  the  shift 
and  spread  observed.  Bottom  left:  SM  cut  for  66  km  without  platform  maneuver. 
Bottom  right:  SM  cut  for  66  km  with  pitch.  The  antenna  pattern  resolution  is  not 
high  enough  to  produce  significant  changes  in  the  PSD  cut.  A  slight  negative  Doppler 
shift  and  shortening  of  the  nnlls  is  observed.  The  nnlls  wonld  eventnally  fill  and  the 
clntter  null  shift  if  more  channels  were  added. 


the  120  X  1  linear  array  and  66km  case,  16  elevation  channels  are  needed  to  reduced 
the  clutter  illumination  at  0  =  0°  3  dB.  That  elevation  channel  addition  increases 
MNP  from  3840  to  61440.  The  slight  redaction  in  clutter  power  is  heavily  trumped 
by  the  enormous  increase  in  required  processing  power,  especially  when  pitch  effects 
are  minimal  at  distant  ranges  for  the  parameters  simulated. 
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Figure  4.11:  SINK  Loss  for  11  km  with  pitch.  The  number  of  azimuth  channels  is 
varied,  while  M  =  20,  P  =  1.  In  the  left  plot,  the  number  of  channels  remains  small 
enough  such  that  the  antenna  pattern  lacks  enough  resolution  for  pitch  effects  to  be 
seen.  As  the  number  of  channels  increases,  the  null  narrows  as  would  be  expected 
if  pitch  were  not  involved.  In  the  right  plot,  the  number  of  channel  is  high  enough 
such  that  pitch  effects  are  observed.  The  null  shifts  and  widens  due  to  the  increased 
azimuth  resolution  and  resulting  PSD  shift  and  spread. 


Azimuth  Angle  Azimuth  Angle 


Figure  4.12:  Antenna  pattern  azimuthal  cut  at  11  km  and  12.86°  pitch.  The 

addition  of  five  (5)  elevation  channels  results  in  a  significant  reduction  in  clutter 
illumination  at  11  km.  The  added  elevation  definition  placed  the  clutter  at  11  km  in 
previously  non-existent  sidelobes.  The  left  plot  is  an  azimuth  cut  for  a  20  x  1  array, 
showing  a  clutter  illumination  at  0  =  0  of  22dB.  The  right  plot  is  an  azimuth  cut  for 
a  20  X  5  array,  showing  a  clutter  illumination  at  0  =  0  of  lOdB. 

Linear  array  pitch  effects  are  very  significant  at  close  ranges.  Pitch  rotation 
changes  the  azimuth  angles  to  clutter  patches  more  quickly  at  close  ranges  than 
distant  ranges.  Multiple  antenna  pattern  lobes  illuminate  each  clutter  patch,  resulting 
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Figure  4.13:  The  addition  of  five  (5)  elevation  channels  places  clutter  at  11km  in 
elevation  sidelobes  not  present  in  the  linear  array.  The  net  result  is  a  significant 
reduction  in  clutter  power.  The  reduction  is  not  observed  at  66km  due  to  the  small 
elevation  angle  to  that  range. 

in  a  broad  Doppler  spread  and  shift.  The  broadness  of  the  shift  is  proportional  to 
the  azimuth  resolution  (and  hence  number  of  azimuth  channels)  in  the  array.  The 
addition  of  elevation  channels  significantly  improves  performance  at  close  ranges  due 
to  the  reduction  of  clutter  illumination  and  constant  target  illnmination.  Compared 
to  the  distant  target  scenario,  fewer  elevation  channels  result  in  significant  clutter 
reduction. 

Pitch  effects  for  the  linear  and  planar  array,  at  11  km  and  66  km,  were  addressed 
in  this  section.  Platform  roll  is  investigated  next. 


92 


Figure  4.14:  MVE  for  20  x  1  and  20  x  5  arrays  with  and  without  pitch  at  11  km. 
The  addition  of  elevation  channels  with  pitch  impacts  the  clutter  spectrum.  The  top 
plots  show  the  MVE  for  a  20  x  1  planar  array  without  (left)  and  with  (right)  pitch  at 
11km.  The  bottom  plots  show  the  MVE  for  the  20  x  5  array  with  (right)  and  without 
(right)  pitch.  Less  spreading  occurs  in  the  20  x  5  case. 

4 -3  Roll  effects 

Platform  roll  corresponds  to  rotation  about  the  x-axis.  The  antenna  is  steered 
to  boresight,  where  the  target  is  located  at  the  start  of  the  CPI.  In  the  roll  scenario, 
the  target  is  not  always  located  at  boresight  due  the  platform  maneuver.  Therefore, 
the  target  is  not  constantly  illuminated  (^^  7^  1m)-  R-oll  effects  will  therefore  be  due 
to  changes  in  both  target  and  clutter  illumination. 


4.3.1  Linear  Array  with  Target  at  11  km  and  66  km.  Eirst  consider  a  target 
at  66km.  Roll  changes  in  the  azimuth  and  elevation  angles  to  the  location  0  =  0 
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Figure  4.15:  Output  SINK  for  20  x  1  and  20  x  5  array  with  pitch.  The  addition  of 
five  (5)  elevation  channels  greatly  improves  SINK  and  narrows  the  clutter  null  to  near 
stationary  MF  performance.  The  left  plot  shows  output  SINK  for  a  20  x  1  array  for 
11km.  The  pitch  causes  a  widening  of  the  clutter  null  around  zero  Doppler.  The  right 
plot  shows  output  SINK  for  a  20  x  5  array  at  11km.  The  null  narrows  and  shortens 
due  to  decreased  clutter  illumination  caused  by  the  elevation  definition  added  and 
added  DOF. 
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Figure  4.16:  A  12.86°  roll  rotation  results  in  no  change  in  azimuth  angle  and  —12.86° 
change  in  elevation  angle  to  the  location  0  =  0  at  a  range  of  66km. 

at  66km  are  depicted  in  Figure  4.16.  Notice  a  12.86°  roll  rotation  corresponds  to 
no  change  in  azimuth  and  an  one-for-one  change  in  elevation  of  —12.86°  for  (p  —  0. 
Significant  SINK  impact  is  therefore  expected  only  if  the  array  pattern’s  elevation 
definition  changes  considerably  over  that  12.86°  span.  In  the  case  of  the  linear  array, 
the  pattern  lacks  elevation  definition.  The  poor  elevation  definition  leads  to  little 
difference  in  clutter  illumination  in  both  the  roll  right  and  roll  left  scenarios,  as 
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Figure  4.17:  Antenna  pattern  azimnthal  cuts  at  66  km  and  12.86°  roll.  On  the 

left,  an  azimuth  cut  for  the  roll  right  scenario  shows  no  difference  than  the  roll  left 
scenario  (right).  The  lack  of  elevation  definition  in  the  linear  array  causes  the  lack  of 
illumination  difference. 

depicted  in  Fignre  4.17.  In  the  case  of  roll  right,  the  mainbeam  is  rolling  skyward 
away  from  the  ground.  In  the  roll  left  scenario,  the  mainbeam  is  rolling  downward 
toward  the  ground.  A  slight  change  of  —0.45  dB  in  target  illumination  occurs  due  to 
roll  in  either  direction.  With  little  change  in  clutter  or  target  illumination  at  66  km, 
roll  has  an  insignificant  impact  on  the  clntter  spectra  and  ontpnt  SINK  as  shown 
Figure  4.18.  Adding  azimuth  channels  narrows  the  clutter  null  and  improves  SINK, 
as  depicted  in  Figure  4.19.  The  clutter  null  narrows  and  SINK  loss  decreases  as  the 
number  of  azimuth  channels  is  increased. 

Snmmarizing  roll  effects  on  the  linear  array,  at  both  11km  and  66km  the  lack 
of  antenna  pattern  elevation  definition  results  in  insignificant  changes  in  target  and 
clutter  illumination.  These  minor  differences  translate  to  minor  differences  in  the 
clntter  spectrnm  and  ontpnt  SINK.  The  addition  of  azimuth  channels  improves  SINK 
and  narrows  the  clntter  nnll. 


4-3.2  Planar  Array  with  Target  at  11  km  and  66  km.  In  the  linear  array, 
the  direction  of  role  was  of  no  consequence  due  to  the  lack  of  elevation  definition  in 
the  antenna  pattern.  The  addition  of  elevation  channels  makes  the  direction  of  roll 
important.  In  the  roll  right  scenario,  the  mainbeam  rolls  skyward  away  from  the 
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Figure  4.18:  MVE,  SNR,  SINK  for  20  x  1  array  at  66  km  with  roll.  The  top  left  hgure 
is  the  MVE  for  the  stationary  platform.  The  top  right  figure  is  the  MVE  with  roll, 
showing  no  discernable  difference  from  the  stationary  case.  The  SNR  is  depicted  in 
the  bottom  left  figure,  illustrating  just  a  —0.45  dB  change  in  target  illumination.  The 
little  change  in  clutter  spectrum  coupled  with  the  small  change  in  target  illumination 
results  in  minimal  change  in  the  output  SINR,  shown  in  the  bottom  right  figure. 


clutter  and  illuminates  the  clutter  with  more  sidelobes.  In  the  roll  left  scenario,  the 
mainbeam  actually  rolls  towards  the  clutter  and  illuminates  more  clutter  with  the 
mainbeam.  The  net  result  is  two  very  different  clutter  spectrums,  depicted  in  Figure 
4.20.  Roll  left  yields  a  14dB  increase  in  clutter  power  because  the  mainbeam  rotates 
toward  the  ground,  while  the  roll  right  yields  a  lOdB  decrease  because  the  mainbeam 
rotates  skyward.  The  target  illumination  also  undergoes  drastic  changes  compared 
to  the  linear  case,  depicted  in  Figure  4.21.  However,  these  dramatic  differences  in 
clutter  spectra  do  not  translate  to  differences  in  output  SINR  depicted  in  Figure  4.22. 
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Figure  4.19:  SINK  loss  and  at  11  km  and  66  km  with  roll.  The  clutter  null  narrows 
and  SINK  loss  decreases  as  the  number  of  azimuth  channel  is  increased  in  both  the 
11  km  (left)  and  66  km  (right)  scenario. 


Despite  a  22dB  difference  in  clutter  spectra,  output  SINK  appears  identical.  The 
result  is  attributable  to  the  high  number  of  MF  DOF,  and  the  fact  that  the  clutter 
spectrum  did  not  shift  in  Doppler. 

Next,  the  relationship  between  increasing  elevation  channels  and  roll  effects  is 
considered.  In  the  linear  array,  adding  azimuth  channels  narrowed  the  clntter  null. 
However,  the  addition  of  elevation  channels  has  the  opposite  effect.  The  SINK  loss 
plot  in  Figure  4.23  shows  a  widening  null  and  increasing  loss  and  the  number  of 
elevation  channels  increase. 

Roll  effects  are  minimal  for  the  linear  array,  at  all  ranges,  due  to  the  lack  of 
elevation  definition  in  the  antenna  pattern.  The  array  remains  in  its  original  place 
throughout  the  CPI,  and  the  change  in  clutter  and  target  illumination  is  minute.  The 
addition  of  elevation  channels  enhances  roll  effects. 


4-4  Yaw  Effects 

Platform  yaw  corresponds  to  rotation  abont  the  z-axis.  The  antenna  is  steered 
to  boresight,  where  the  target  is  located  at  the  start  of  the  CPI.  In  the  yaw  scenario, 
the  target  is  not  always  located  at  boresight  due  the  platform  maneuver.  Therefore, 
the  target  is  not  constantly  illnminated  (^^  7^  1m)- 
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Figure  4.20:  MVE  and  SM  cut  for  planar  arrays  with  and  without  roll.  The  addition 
of  elevation  channels  with  roll  impacts  the  clutter  spectrum  at  both  11km  and  66km. 
The  top  plots  show  the  MVE  for  a  20  x  5  planar  array  with  roll  left  (left)  and  roll 
right  (right)  at  11km.  A  broadening  in  Doppler  is  observed  in  both  cases,  but  no 
magnitude  information  can  be  taken  from  MVE.  The  bottom  plot  shows  the  SM  PSD 
cuts  at  zero  azimuth  for  roll  right,  roll  left,  and  stationary.  Roll  left  yields  a  14dB 
increase  in  clutter  power  because  the  mainbeam  rotates  toward  the  ground,  while  the 
roll  right  yields  a  lOdB  decrease. 

Jt-.Ji-.l  Linear  Array  with  Target  at  11  km  and  66  km.  Yaw  changes  in  the 
azimuth  and  elevation  angles  to  the  location  0  =  0  at  66km  are  depicted  in  Eigure 
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Figure  4.21:  Signal-to-Noise  ratio  for  linear  and  planar  array  and  target  at  11km 

with  roll.  The  poor  elevation  definition  of  the  20  x  1  linear  array  (left)  yields  little 
change  in  target  illumination  despite  12.86°  of  roll.  The  addition  of  five  (5)  eleva¬ 
tion  channels  (right)  increases  elevation  definition  and  yields  a  — 7dB  swing  in  target 
illumination. 


Figure  4.22:  SINK  for  20  X  5  planar  array  with  roll  at  11  km.  Left:  roll  right. 

Right:  roll  left.  The  increased  clutter  power  is  evident  in  the  roll  left  plot  in  the  large 
drop  in  NA  performance.  However,  the  MF  overcomes  the  increase  in  clutter  with 
MNP  =  3200  DOF. 

4.24.  Notice  a  12.86°  yaw  rotation  corresponds  to  no  change  in  elevation  and  an 
one-for-one  change  in  azimuth  of  —12.86°  for  0  =  0.  This  relationship  is  the  same  for 
11km  and  66km.  Significant  SINK  impact  is  therefore  expected  if  the  array  pattern’s 
azimuth  definition  changes  considerably  over  that  12.86°  span.  In  the  case  of  the  linear 
array,  the  pattern  has  significant  azimuth  definition  for  any  realistic  array  length.  The 
azimuth  definition  leads  to  significant  differences  in  the  clutter  illumination  for  both 


99 


Figure  4.23:  SINK  loss  varying  elevation  channels  with  roll.  The  left  plot  shows 

SINK  loss  for  the  11  km  scenario  and  increasing  elevation  channels  and  the  left  plot 
shows  the  same  for  66  km.  In  both  cases,  the  loss  and  null  width  increase  as  P 
increases.  N  =  20  for  all  plots. 


Figure  4.24:  A  12.86°  yaw  rotation  results  in  no  change  in  elevation  angle  and 

—  12.86°  change  in  azimuth  angle  to  the  location  0  =  0  at  both  11km  and  66km. 

11km  and  66km,  as  shown  in  Figure  4.25.  This  same  azimuth  definition  leads  to 
significant  differences  in  the  clutter  spectrum  for  both  11km  and  66km,  as  shown 
in  Figure  4.26.  Smearing  in  both  azimuth  and  Doppler  is  apparent  in  Figures  4.26. 
However,  the  SM  PSD  cut  does  not  look  drastically  changed  from  the  stationary 
case.  As  the  platform  yaws  clockwise  (CW)  about  the  z-axis,  the  mainbeam  and 
antenna  pattern  rotates  as  well.  The  velocity  vector  and  consequently  the  Doppler 
of  each  patch  changes.  However,  the  overall  magnitude  of  Doppler  illumination 
remains  relatively  unchanged.  In  the  stationary  case,  the  patch  at  0  =  0  had  a 
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Figure  4.25:  Antenna  pattern  azimuthal  cuts  for  20  x  1  at  11  and  66  km  with 

yaw.  Yaw  rotation  about  the  z-axis  produces  a  one-for-one  change  in  azimuth  angle. 
A  12.86°  rotation  corresponds  to  nearly  two  lobes  (mainbeam  plus  sidelobe)  of  the 
pattern. 


Doppler  frequency  of  zero  and  was  at  boresight  throughout  the  CPI.  Now  assume  the 
platform  yawed  1°.  That  same  patch  is  no  longer  at  at  boresight  nor  is  its  Doppler 
frequency  zero.  However,  the  patch  previously  at  0  =  1°  is  now  at  boresight  and  its 
Doppler  frequency  is  zero.  When  the  platforms  yaws  2°,  the  patch  previously  at  2° 
and  non-zero  Doppler  is  now  at  boresight  and  zero  Doppler.  The  radar  is  “looking”  at 
its  boresight  throughout  the  CPI,  and  that  view  is  not  drastically  changing.  However, 
a  yaw  has  an  enormous  impact  on  target  illumination,  depicted  in  Figure  4.27.  The 
12.86°  yaw  puts  the  target  though  the  hrst  antenna  pattern  null  and  halfway  into 
the  second  null.  The  target  is  nearly  begin  “yawed”  out  of  the  picture.  The  result 
is  a  significant  degradation  in  output  SINK,  shown  in  Figure  4.28  for  both  11  km 
and  66  km.  If  the  number  of  azimuth  channels  is  increased,  the  antenna  pattern 
beamwidth  decreases  and  the  target  illumination  change  becomes  more  pronounced. 
The  SINK  loss  plot  in  Figure  4.29  shows  loss  proportional  to  the  number  of  azimuth 
channels  N .  The  larger  N  becomes,  the  close  the  radar  gets  to  “yawing”  the  target 
out  of  the  picture. 


4.4--2  Planar  Array  with  Target  at  11  km  and  66  km.  The  addition  of 
elevation  elements  once  again  changes  the  way  the  clutter  is  illuminated.  The  resulting 
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Figure  4.26:  Clutter  spectrum  for  20  x  1  arrays  with  and  without  yaw  at  11  km 

and  66  km.  The  top  plots  show  the  MVE  for  11  km  (left)  and  66  km  (right)  with 
yaw.  These  plots  show  extreme  spread  in  Doppler  and  azimuth.  The  bottom  plots 
show  the  SM  PSD  cut  at  zero  azimuth  for  11  km  (left)  and  66  km  (right).  Unlike  the 
MVE,  drastic  differences  are  not  seen  in  the  SM  cut. 

clutter  spectrum  are  depicted  for  the  20  x  5  array  at  11  km  and  66  km  in  Eigure 
4.30.  Target  illumination  is  shown  in  Eigure  4.31.  The  addition  of  elevation  channels 
does  not  impact  the  target  illumination,  as  the  target  illumination  depends  only  on 
the  azimuthal  definition  due  to  yaw  being  a  rotation  about  the  z-axis.  However, 
the  target  illumination  remains  drastically  reduced  compared  to  the  stationary  case. 
As  the  number  of  elevation  channels  is  increased,  the  added  DOE  and  reduction  in 
clutter  illumination  do  not  overcome  the  drastic  reduction  in  target  illumination. 
Some  improvement  is  observed  however  at  the  11  km  range,  as  the  change  in  clutter 
illumination  impacts  performance  more  so  than  at  66  km.  Eigure  4.32  shows  the 
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Figure  4.27:  SNR  for  20  x  1  array  with  yaw  at  11  km  and  66  km.  The  yaw  induces 
major  changes  in  target  illumination  in  both  the  11  km  (left)  and  66  km  (right) 
scenarios.  The  target  is  illuminated  nearly  6  dB  less  only  a  third  of  the  way  through 
the  CPI.  By  the  halfway  point  of  the  CPI,  the  difference  is  30  dB  down.  The  12.86° 
yaw  puts  the  target  though  the  first  antenna  pattern  null  and  halfway  into  the  second 
null. 


Figure  4.28:  Output  SINK  for  20  x  1  array  with  pitch  at  11  km  and  66  km.  The 
vast  reduction  in  target  illumination  results  in  significant  impact  on  SINK  for  both 
the  11  km  (left)  and  66  km  (right)  scenario. 

SINK  loss  for  varying  elevation  channels.  In  the  11  km  scenario,  a  slight  narrowing 
of  the  clutter  null  is  observed  and  loss  is  inversely  proportional  to  P.  It  is  noted 
however  that  the  SINK  loss  maximum  is  -12dB,  despite  the  narrowing  null.  While 
the  decrease  in  clutter  illumination  lessens  the  degradation  caused  by  yaw,  the  hit  it 
target  illumination  can  not  be  overcome.  In  the  66  km  case,  adding  elevation  channels 
does  not  effect  the  null  width  or  loss.  Physically,  the  66  km  range  is  farther  away 
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Figure  4.29:  SINK  loss  at  11  and  66  km  with  roll.  Severe  loss  is  observed  at  both 
range,  11  km  (left)  and  66  km  (right).  The  beamwidth  associated  with  these  values 
of  N  is  small  enough  such  that  the  target  illumination  is  drastically  reduced. 

where  the  added  elevation  channels  do  not  effect  the  clutter  illumination.  The  result 
is  again  limited  by  the  decreased  target  illumination. 

Overall,  yaw  is  by  far  the  most  devastating  maneuver  effect,  at  short  and  long 
ranges.  The  main  culprit  is  drastically  lower  power  with  which  the  target  is  illu¬ 
minated.  The  addition  of  elevation  channel  is  ineffective  at  mitigating  the  impact. 
Pitch  contributions  are  limited  unless  the  number  of  azimuth  channels  produce  a 
beamwidth  on  the  order  of  the  change  in  azimuth  and  elevation  to  a  range.  The 
number  of  azimuth  channels  needed  is  proportional  to  the  range  of  interest.  The 
addition  of  elevation  channel  mitigated  pitch  effects.  Roll  effects  were  negligible  for 
the  linear  array  due  to  lack  of  elevation  definition.  Adding  azimuth  channels  to  the 
linear  array  improved  performance.  The  addition  of  elevation  channels  enhanced  roll 
effects,  as  the  added  elevation  definition  resulted  in  fluctuations  in  target  and  clutter 
illumination  not  seen  with  the  linear  array. 

4-5  Maneuver  Effeets  on  Sample  Support  Homogeneity 

The  statistical  clutter  environment  is  rarely,  if  ever,  known  a  priori  and  must 
be  estimated  based  on  snapshots  from  range  cells  other  than  the  RUT.  An  accurate 
estimate  relies  upon  homogeneous,  independent-identically  distributed  (iid)  sample 


104 


Figure  4.30:  MVE  for  20  x  5  arrays  with  and  without  yaw  at  11  km  and  66  km.  The 
MVE  for  11  km  in  the  top  right  plot  shows  severe  smearing  in  azimuth  and  Doppler. 
Similar  smearing  is  seen  in  the  MVE  at  66  km  (top  right).  The  SM  cuts  at  11  km 
(bottom  right)  and  66  km  (bottom  left)  show  a  shift  towards  negative  Doppler. 

support  vectors.  For  this  research,  the  generalized  inner-product  [11], 


GIPk  —  xfR-RUTXk- 


(4.15) 


is  used  to  access  maneuver  induced  heterogeneities  in  sample  support  data.  Using 
the  known  covariance  of  the  RUT,  the  GIF  compares  the  magnitude  and  phase  of 
snapshots  at  surronnding  range  cells.  The  larger  the  GIF  variance  and  range  across 
a  set  of  sample  support  vectors,  the  greater  the  non- homogeneity.  The  key  physical 
insight  into  what  the  GIF  does  is  that  it  measures  the  statistical  “likeness”  of  the  RUT 
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Figure  4.31:  SNR  for  20  x  5  array  with  yaw  at  11  km  and  66  km.  The  addition  of 
elevation  does  not  impact  the  major  changes  in  target  illumination  in  both  the  11  km 
(left)  and  66  km  (right)  scenarios.  The  target  is  illumination  nearly  6  dB  down  only 
a  third  of  the  way  through  the  CPI.  By  the  halfway  point  of  the  CPI,  the  difference 
is  30  dB  down.  The  12.86°  yaw  puts  the  target  though  the  first  antenna  pattern  null 
and  halfway  into  the  second  null. 


Figure  4.32:  SINK  loss  varying  elevation  channels  at  11  and  66  km  with  yaw.  At 

11  km  (left),  a  slight  narrowing  of  the  clutter  null  is  observed  and  loss  is  inversely 
proportional  to  P.  The  SINK  loss  maximum  is  -12dB,  despite  the  narrowing  null. 
While  the  decrease  in  clutter  illumination  lessens  the  degradation  caused  by  yaw 
at  11  km,  the  hit  it  target  illumination  can  not  be  overcome.  At  66  km  (right), 
adding  elevation  channels  does  not  effect  the  null  width  or  loss.  Physically,  the  66  km 
range  is  farther  away  where  the  added  elevation  channels  does  not  effect  the  clutter 
illumination.  The  result  is  again  limited  by  the  decreased  target  illumination. 

with  other  range  cells.  The  more  statistically  alike  they  are,  the  more  homogenous 
that  sample  support  vector  is. 
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Figure  4.33:  GIF  for  20  X  1,  known  covariance  at  66  km  and  no  maneuver.  The 

CPI  length  is  20  pulses.  Note  the  value  at  the  RUT  (550)  is  26  dB,  corresponding  to 
the  dB  value  of  MNP.  A  total  of  338  range  cells  are  within  3dB  of  the  RUT. 

Sample  support  homogeneity  is  examined  for  a  20  x  1  linear  and  20  x  4  planar 
array.  Recalling  the  sample  support  discussion  from  Chapter  II,  Factored-Time  Space 
(FTS)  requires  2N  i.i.d.  sample  support  vectors  and  Joint-Domain-Localized  (JDL) 
requires  2riaribrie  to  achieve  within  3dB  of  known  covariance  [2] .  FTS  therefore  requires 
40  i.i.d.  sample  support  vectors.  Note,  the  terminology  “3dB  of  known  covariance” 
means  known  FTS  and  JDL  covariance.  JDL  may  be  performed  with  fewer  DOF,  but 
is  not  assigned  a  minimum  requirement  in  this  exercise.  As  the  maneuver  induced 
heterogeneities  are  accessed,  these  i.i.d.  sample  support  requirements  are  kept  in 
mind.  First,  the  GIP  results  for  the  homogeneous  cases  are  examined.  The  CPI  is 
shorted  to  20  pulses  for  all  GIP  simulations  in  the  interest  of  computational  loads. 
In  Figure  4.34,  the  RUT  is  range  cell  550  at  66  km.  The  GIP  value  at  the  RUT 
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is  26dB,  corresponding  to  the  dB  value  of  MNP.  In  general,  the  farther  away  a 
range  cell  is  from  the  RUT,  the  less  statistically  alike  the  two  are  and  the  larger 
the  GIP.  From  a  physical  standpoint,  the  range  cells  at  close  ranges  are  under  much 
higher  illumination  and  the  phase  of  snapshots  at  close  ranges  differ  greatly  from  the 
RUT.  Figure  4.34  illustrates  this  concept  nicely.  No  published  work  exists  stating  a 
quantitative  maximum  GIP  that  equates  to  the  sample  support  vectors  under  that 
maximum  are  “homogeneous  enough”  to  meet  the  i.i.d.  assumption.  The  focus  of 
this  research  is  not  determining  which  range  cells  are  intelligent  choices  for  sample 
support.  Instead,  the  focus  is  what  effects  does  platform  maneuver  have  on  the  sample 
support  homogeneity.  For  comparisons  between  range  cells,  a  3dB  threshold  is  used 
in  this  research.  In  the  stationary  case  of  Figure  4.34,  a  total  of  338  range  cells  (from 
number  292  and  higher)  are  within  3dB  of  the  RUT.  Sufficient  homogeneous  sample 
support  is  therefore  available  for  both  FTS  and  JDL. 

Range  cell  292  is  30.9km  from  the  RUT.  The  likelihood  of  range  cells  of  a 
statistically  “alike”  nature  existing  at  ranges  farther  away  than  that  is  extremely 
low.  The  swath  of  range  cells  investigated  is  therefore  limited  to  number  292  and 
higher.  The  addition  of  platform  maneuver  is  shown  in  Figure  ??.  The  addition  of 
roll  had  minimal  effect  on  the  GIP.  All  338  range  cells  within  3dB  of  the  RUT  remain 
under  3dB.  However,  pitch  reduced  that  number  to  271  and  yaw  reduced  it  to  278. 
Overall,  maneuver  effects  on  sample  support  homogeneity  is  not  dramatic.  The  worst 
case  (pitch)  yields  a  20%  reduction  of  range  cells  within  3dB,  but  the  total  number 
of  sample  support  vectors  within  3dB  remains  high  enough  for  partially  adaptive 
processing  requirement. 

The  stationary  GIP  with  the  addition  of  four  (4)  elevation  channels  is  shown  in 
Figure  4.35.  Note  the  reduction  in  the  number  of  sample  support  vectors  within  3dB 
when  elevation  channels  are  added.  There  are  now  258  range  cells  within  3dB  of  the 
RUT.  The  addition  of  elevation  channel  means  added  definition  to  the  antenna  pat¬ 
tern.  Therefore,  the  differences  in  illumination  between  any  two  range  cells  increases. 
The  addition  of  platform  maneuver  is  shown  in  Figure  4.36.  Overall,  the  maneuver 
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Figure  4.34:  GIF  for  20  x  1,  known  covariance  at  66  km  and  with.  The  CPI  length 
is  20  pulses.  In  the  stationary  and  roll  cases,  338  rangge  cells  are  within  3dB  of  the 
RUT  GIP.  However,  pitch  has  reduced  that  number  to  271  and  yaw  to  278. 

had  little  impact  on  the  number  of  sample  support  vectors  within  3dB  of  the  RUT. 
Pitch  yielded  253  (98%  of  stationary  total),  roll  yielded  264  (102%),  and  yaw  yielded 
247  (96%).  Overall  for  a  target  at  66  km,  platform  maneuver  has  little  impact  on 
data  homogeneity  around  the  RUT.  The  changes  in  angles  at  distant  ranges  is  small, 
as  are  the  resulting  changes  in  illumination  and  phase  of  the  returns. 

Attention  is  now  turned  to  a  target  at  11  km.  Unlike  the  66  km  scenario, 
platform  maneuver  introduces  significant  changes  in  angle  and  illumination  from  one 
location  to  another,  as  shown  in  Section  4. 2-4. 4.  The  stationary  GIP  is  shown  in 
Figure  4.37.  Note  the  value  at  the  RUT  (92)  is  26  dB,  corresponding  to  the  dB  value 
of  MNP.  A  total  of  592  range  cells  are  within  3dB  of  the  RUT.  The  maneuver  results 
are  shown  in  Figure  4.38.  Pitch  yielded  450  (76%  of  stationary  total),  roll  yielded  596 
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Figure  4.35:  GIF  for  20  X  4,  known  covariance  at  66  km  and  no  maneuver.  The 

CPI  length  is  20  pulses.  Note  the  value  at  the  RUT  (550)  is  32  dB,  corresponding  to 
the  dB  value  of  MNP.  A  total  of  258  range  cells  are  within  3dB  of  the  RUT. 

(96%),  and  yaw  yielded  393  (66%)  range  cells  within  3dB  of  the  RUT.  However,  range 
cells  around  the  RUT  are  highly  corrupted.  Covariance  estimation  is  often  performed 
with  sample  support  vectors  near  ranges  of  the  RUT.  However,  maneuver  makes  these 
surrounding  range  cells  highly  heterogenous. 

Overall  for  the  linear  array,  sample  support  homogeneity  effects  are  increased  at 
closer  range.  At  66  km,  sample  support  within  3dB  of  the  RUT  is  80%  of  the  station¬ 
ary  total  when  pitch  is  introduced,  100%  for  roll,  and  82%  for  yaw.  At  11  km,  those 
numbers  fall  to  76%,  96%,  and  66%  respectively.  Additionally,  the  cells  surrounding 
the  RUT  are  highly  heterogenous  at  11  km.  Sample  support  homogeneity  increases 
when  elevation  channels  are  added,  for  both  11  km  and  66  km.  The  increased  el¬ 
evation  channels  decreases  the  overall  illumination  of  the  clutter.  When  maneuver 
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Figure  4.36:  GIF  for  20  X  4,  known  covariance  at  66  km  with  maneuver.  The  CPI 
length  is  20  pulses  with  the  RUT  =550.  Pitch  yielded  253  (98%  of  stationary  total), 
roll  yielded  264  (102%),  and  yaw  yielded  247  (96%)  range  cells  within  3dB  of  the 
RUT. 

is  introduced,  the  changes  is  illumination  are  lessened  and  the  resulting  data  more 
“alike”. 

4.-6  Results  Summary 

Pitch  maneuver  effects  are  limited  unless  the  number  of  azimuth  channels  pro¬ 
duce  a  beamwidth  on  the  order  of  the  change  in  azimuth  and  elevation  to  the  range 
of  interest.  MF  performance  improves  as  the  number  of  azimuth  channels  increases, 
until  that  threshold  pattern  resolution  is  achieved.  The  number  of  azimuth  channels 
needed  is  proportional  to  the  range  of  interest,  thus  pitch  effects  are  strongest  at 
close  ranges.  The  addition  of  elevation  channel  mitigate  pitch  effects,  but  the  added 
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Figure  4.37:  GIF  for  20  X  1,  known  covariance  at  11  km  and  no  maneuver.  The 

CPI  length  is  20  pulses.  Note  the  value  at  the  RUT  (92)  is  26  dB,  corresponding  to 
the  dB  value  of  MNP.  A  total  of  592  range  cells  are  within  3dB  of  the  RUT. 

DOF  are  unnecessary  at  ranges  where  pitch  is  a  non-factor.  The  target  resides  in  the 
mainbeam  at  Foresight  throughout  the  CPI,  due  to  the  rotation  about  the  Foresight 
axis. 

Roll  has  little  impact  on  the  linear  array  at  all  ranges,  due  the  the  lack  of 
elevation  definition  in  the  antenna  pattern.  That  lack  of  definition  translates  to  min¬ 
imal  changes  in  clutter  and  target  illumination.  Increases  in  the  number  of  azimuth 
channels  narrows  the  clutter  null  and  improves  MF  performance.  The  addition  of 
elevation  channels,  however,  magnifies  roll  impact.  The  added  elevation  definition 
changes  clutter  and  target  illumination,  resulting  in  degraded  performance. 
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Figure  4.38:  GIF  for  20  X  1,  known  covariance  at  11  km  with  maneuver.  The  CPI 
length  is  20  pulses  with  the  RUT  =92.  Pitch  yielded  450  (76%  of  stationary  total), 
roll  yielded  596  (96%),  and  yaw  yielded  393  (66%)  range  cells  within  3dB  of  the  RUT. 
Range  cells  around  the  RUT  are  very  heterogenous. 

The  most  devastating  maneuver  effect  at  all  ranges  is  yaw.  The  main  culprit  is 
drastically  lower  power  with  which  the  target  is  illuminated.  Increasing  the  number 
of  azimuth  channels  tightens  the  antenna  beamwidth,  leading  to  decreased  target  illu¬ 
mination  and  MF  performance.  The  addition  of  elevation  channels  slightly  improves 
performance  at  close  ranges,  where  the  decrease  in  clutter  illumination  is  noticeable. 
As  range  increases,  the  added  elevation  channels  have  less  impact  as  the  clutter  illu¬ 
mination  decreases  with  range. 

Similar  to  the  MF  SINR  results,  platform  maneuver  impacts  on  sample  support 
homogeneity  are  more  pronounced  at  11km  than  66km.  The  changes  in  azimuth 
and  elevation,  and  therefore  the  changes  in  the  magnitude  and  phase  of  the  clutter 
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Table  4.2:  Percent  reduction  in  number  of  sample 

support  vectors  GIP  falling  within  3dB  of  the  RUT. 


Pitch 

Roll 

Yaw 

11  km 

66  km 

11  km 

66  km 

11  km 

66  km 

Linear 

24% 

20% 

4% 

0% 

34% 

18% 

Planar 

4% 

2% 

3% 

0% 

6% 

4% 

snapshots,  are  increased.  The  largest  impact  is  observed  in  the  linear  array  with  yaw 
at  11km,  where  only  66%  of  sample  support  vectors  originally  within  3dB  of  the  RUT 
remain  within  3dB.  This,  and  all  GIP  results,  are  highly  optimistic  given  the  array 
sizes  tested.  Recall  pitch  effects  are  not  seen  in  the  output  SINR  plots  for  66km  until 
N  approached  115  elements.  Yet  20%  of  sample  support  previously  with  3dB  of  the 
RUT  was  lost  when  the  array  was  just  20  x  1. 
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V.  Conclusions 


This  research  develops  a  space-time-adaptive-processing  (STAP)  model  incorpo¬ 
rating  platform  pitch,  yaw,  and  roll.  It  characterizes  the  resulting  interference 
environment,  and  examines  the  maneuver  impact  on  fully  adaptive  processing  perfor¬ 
mance  and  homogeneous  sample  support.  Chapter  II  summarizes  the  airborne  radar 
problem,  the  stationary  STAP  model  of  [3,10],  and  STAP  fundamentals.  The  model 
of  [3, 10]  serves  as  the  foundation  for  the  maneuver  model.  Chapter  III  details  the 
maneuver  model,  including  modified  signal,  target,  clutter,  and  formatting  models. 
Chapter  IV  simulates  Chapter  Ill’s  model.  Each  maneuver  and  its  effect  on  the  in¬ 
terference  environment,  target,  Signal-to-Interference-plus-Noise- Ration  (Ratio),  and 
homogeneous  sample  support  are  examined  in  detail. 

5.1  STAP  Maneuver  Model  for  Side-Looking  Arrays 

The  maneuver  model  developed  can  simulate  any  degree  of  pitch,  roll,  yaw,  or 
any  combination  of  the  three.  This  research  focused  on  development  of  the  model 
itself,  versus  specific  innovative  and  complex  scenarios  to  simulate.  The  framework 
developed,  however,  is  capable  of  handling  complex  rotations  not  addressed  in  this 
research’s  simulations.  The  overriding  theme  of  the  maneuver  model  is  that  singular 
quantities  in  the  stationary  model  of  [3,10]  possess  M  unique  quantities  in  the  maneu¬ 
ver  model.  The  platform  defines  the  reference  coordinate  system  at  all  times.  The  {9, 
(j))  radar  coordinates  of  each  clutter  patch  and  the  target  are  first  transformed  into 
their  cartesian  counterparts  (x,y,z).  Rotation  matrices  are  used  to  recalculate  the 
coordinates  to  each  {x,y,z)  point  based  on  the  amount  of  rotation,  measured  with  the 
platform  orientation  as  the  reference.  The  resulting  {x' ,y' ,z')  coordinates  are  trans¬ 
formed  back  to  radar  coordinates  {O' ,  <p').  Each  location  has  M  {9',  (j)')  pairs,  trans¬ 
lating  to  M  Doppler  frequencies,  spatial  frequencies,  and  amplitudes.  Four  MNP  x  1 
vectors  are  created,  one  for  amplitude,  one  each  for  phase  contributions  due  horizontal 
spacing,  vertical  spacing,  and  Doppler.  The  Kronecker  product  paradigm  is  replaced 
with  the  Hadamard  product  in  creating  the  space-time  snapshots. 
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5.2  Impact  on  Matched  Filter  Performance 

Pitch  maneuver  effects  are  limited  unless  the  number  of  azimuth  channels  pro¬ 
duce  a  beamwidth  on  the  order  of  the  change  in  azimuth  and  elevation  to  the  range 
of  interest.  MF  performance  improves  as  the  number  of  azimuth  channels  increases, 
until  that  threshold  pattern  resolution  is  achieved.  The  number  of  azimuth  channels 
needed  is  proportional  to  the  range  of  interest,  thus  pitch  effects  are  strongest  at 
close  ranges.  The  addition  of  elevation  channels  mitigates  pitch  effects,  but  the  added 
DOF  are  unnecessary  at  ranges  where  pitch  is  a  non-factor.  The  target  resides  in  the 
mainbeam  at  boresight  throughout  the  CPI,  due  to  the  rotation  about  the  boresight 
axis. 

Roll  has  little  impact  on  the  linear  array  at  all  ranges,  due  the  the  lack  of 
elevation  definition  in  the  antenna  pattern.  That  lack  of  definition  translates  to  min¬ 
imal  changes  in  clutter  and  target  illumination.  Increases  in  the  number  of  azimuth 
channels  narrows  the  clutter  null  and  improves  MF  performance.  The  addition  of 
elevation  channels,  however,  magnifies  roll  impact.  The  added  elevation  definition 
changes  clutter  and  target  illumination,  resulting  in  degraded  performance. 

The  most  devastating  maneuver  effect  at  all  ranges  is  yaw.  The  main  culprit  is 
drastically  lower  power  with  which  the  target  is  illuminated.  Increasing  the  number 
of  azimuth  channels  tightens  the  antenna  beamwidth,  leading  to  decreased  target  illu¬ 
mination  and  MF  performance.  The  addition  of  elevation  channels  slightly  improves 
performance  at  close  ranges,  where  the  decrease  in  clutter  illumination  is  noticeable. 
As  range  increases,  the  added  elevation  channels  have  less  impact  as  the  clutter  illu¬ 
mination  decreases  with  range. 

5.3  Impact  on  Sample  Support  Homogeneity 

Similar  to  the  MF  SINK  results,  platform  maneuver  impacts  are  more  pro¬ 
nounced  at  11km  than  66km.  The  changes  in  azimuth  and  elevation,  and  therefore 
the  changes  in  the  magnitude  and  phase  of  the  clutter  snapshots,  are  increased.  The 
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largest  impact  is  observed  in  the  linear  array  with  yaw  at  11km,  where  34%  of  sample 
support  vectors  originally  within  3dB  of  the  RUT  are  lost.  This,  and  all  GIF  results, 
are  highly  optimistic  given  the  array  sizes  tested.  Recall  pitch  effects  are  not  seen  in 
the  output  SINR  plots  for  66km  until  N  approached  115  elements.  Yet  20%  of  sample 
support  previously  with  3dB  of  the  RUT  was  lost  when  the  array  was  just  20  x  1. 


5.4  Improvements  and  Future  Research 

The  model  developed  offers  many  potential  extensions  for  future  research.  The 
first  improvement  is  the  addition  of  motion  compensated  steering.  Motion  compen¬ 
sated  steering  contains  two  parts:  compensation  on  transmit  and  compensation  on 
receive.  On  transmit,  simulations  for  this  research  restricted  antenna  steering  to  its 
boresight,  effectively  making  the  radar  unaware  of  its  own  maneuver.  Radar  platforms 
today  have  the  ability  to  monitor  their  position  and  steer  their  transmit  beam  to  a 
fixed  location.  Adding  this  feature  to  the  model  requires  a  simple  modification.  The 
antenna  pattern  presented  as  steered  to  (6*,0)  =  (0,0)  in  Equation  (2.30)  as, 

P-l  JV-l 

iu(d,0)  =  EE  ^j2TTn'dx 

p={)  n=0 


(5.1) 


Steering  the  beam  on  transmit  changes  this  expression  to 


P-l  N-\ 

w[»a)  =  Y.Y. 


^j2TTn'dx  ^j27Tpi^z  ^j27Tn'&l  ^j2TTp^\ 


(6.2) 


where  and  are  the  spatial  frequencies  corresponding  to  look  location  {6\  cjJ'). 
The  target  location  is  therefore  calculated  based  on  the  original  stationary  geometry, 
and  the  M  values  to  that  target  are  calculated  and  create  M  x  1  vectors  of  antenna 
steer.  Ideally,  the  radar  is  able  to  update  its  look  direction  every  range  cell.  The 
second  part  of  motion  compensated  steering,  compensation  on  receive,  is  already 
completely  within  the  context  of  the  model.  Motion  compensation  on  receive  refers 
to  adapting  with  the  motion  compensated  steering  vector.  This  model  calculates 
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that  steering  vector,  but  adapted  with  the  original  uncompensated  steering  vector. 
Mathematically,  the  MF  weight  vector  calculation  changes  from, 

w  =  R“^v,  (5.3) 

to 

w  =  R“^v,  (5.4) 

for  motion  compensated  steering  on  receive. 

The  second  improvement  and  future  extension  of  this  research  is  conversion  from 
SL  to  forward-looking  (FL).  The  conversion  itself  is  rather  simple.  Converting  this 
model  to  FL  requires  a  new  definition  of  the  aircraft  velocity  vector.  It  was  previously 
derived  as. 


Va  =  U^X. 


(5.5) 


The  FL  velocity  vector  is 


Va  =  VaJ. 


(5.6) 


And  the  normalized  Doppler  frequency  of  a  clutter  patch  at  changes  from 


oj  —  2 


Va  cos  9  sin  cj) 


(5.7) 


to 


uj  —  2 


Va  COS  9  COS  (j) 
A/r 


(5.8) 


Additional  modification  is  necessary  in  the  target  model,  as  the  projection  of  the 
aircraft  velocity  vector  to  the  target  location  changes  as  well.  Nonetheless,  the  im¬ 
provement  itself  is  simple  while  the  potential  applications  are  not. 
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The  final  suggestion  for  improvement  and  research  continuation  is  the  use  of 
sub- arrays.  Real  world  antenna  arrays  often  number  in  the  hundreds  or  thousands  of 
elements.  Simulations  for  this  research  were  limited  in  array  size  due  to  computational 
limits.  The  use  of  sub-arrays  removes  that  limitation  by  simulating  the  increased 
array  size  on  transmit,  while  keeping  the  temporal  DOF  at  a  manageable  level.  The 
added  capability  of  sub-arrays  enhances  the  potential  applications  of  this  research 
immensely. 
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